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rvj Abstract 

?H We solve a spectral and an inverse spectral problem arising in the 

V^ computation of peakon solutions to the two-component PDE derived by 

Geng and Xue as a generalization of the Novikov and Degasperis-Procesi 
--. equations. Like the spectral problems for those equations, this one is of a 

'discrete cubic string' type - a nonselfadjoint generalization of a classical 
I I inhomogeneous string - but presents some interesting novel features: there 

Qi^ are two Lax pairs, both of which contribute to the correct complete spectral 

^0 data, and the solution to the inverse problem can be expressed using 

quantities related to Cauchy biorthogonal polynomials with two different 
spectral measures. The latter extends the range of previous applications 
-, of Cauchy biorthogonal polynomials to peakons, which featured either two 

r^ identical, or two closely related, measures. The method used to solve the 

' ' spectral problem hinges on the hidden presence of oscillatory kernels of 

Gantmacher-Krein type implying that the spectrum of the boundary value 
►^ problem is positive and simple. The inverse spectral problem is solved 

_^ by a method which generalizes, to a nonselfadjoint case, M. G. Krein's 

|y-\ solution of the inverse problem for the Stieltjes string. 

00 
/ 1 Introduction 

^^ In this paper, we solve an inverse spectral problem which appears in the context 

^-H of computing explicit solutions to a two-component integrable PDE in 1 + 1 

dimensions found by Geng and Xue [14]. Denoting the two unknown functions 
by u(x,t) and w(a;,i), and introducing the auxihary quantities 






> 

X 



m^ u-Uxx, n^v~Vxx, (1-1) 

we can write the Geng-Xue equation as 



mt + {vrLxU + 3mUx)v — 0, 
TT't + {i^xV + SnVx)u = 0. 



(1.2) 
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(Subscripts denote partial derivatives, as usual.) This system arises as the 
compatibility condition of a Lax pair with spectral parameter z, 



(1.3a) 






vumz ^"2 , (1.3b) 



but also (because of the symmetry in (1.2)) as the compatibility condition of a 
different Lax pair obtained by interchanging u and v, 



(1.4a) 



(1.4b) 

The subject of our paper is the inverse problem of recovering m and n from 
spectral data obtained by imposing suitable boundary conditions on equations 
(1.3a) and (1.4a), in the case when m and n are both discrete measures (finite 
linear combinations of Dirac deltas) with disjoint supports. To explain why this 
is of interest, we will give a short historical background. 

When u = V (and consequently also m = n), the Lax pairs above reduce to 
the Lax pair found by Hone and Wang for V. Novikov's integrable PDE [28, 18] 

rrit + rrixU +3muUx — 0, (1-5) 

and it was by generalizing that Lax pair to (1.3) that Geng and Xue came up 
with their new integrable PDE (1.2). Novikov's equation, in turn, was found as 
a cubically nonlinear counterpart to some previously known integrable PDEs 
with quadratic nonlinearities, namely the Camassa-Holm equation [8] 

fTT-t + rrixU + 2mUx — (1-6) 

and the Degasperis-Procesi equation [10, 9] 

'"Tl't + TTlj-U + imUx = 0. (1-7) 

The equations (1.6) and (1.7) have been much studied in the literature, and 
the references are far too numerous to survey here. Novikov's equation (1.5) is 
also beginning to attract attention; see [16, 17, 19, 21, 25, 27, 30]. What these 
equations have in common is that they admit weak solutions called peakons 
(peaked solitons), taking the form 

N 

w(x,i) = ^77ifc(t)e-l^-^''(*)l, (1.8) 

where the functions Xk{t) and mk{t) satisfy an integrable system of 2 A'' ODEs, 
whose general solution can be written down explicitly in terms of elementary 



functions with the help of inverse spectral techniques. In the Camassa-Holm 

case, this involves very classical mathematics surrounding the inverse spectral 
theory of the vibrating string with mass density g{y), whose eigenmodes are 
determined by the Dirichlet problem 

-(t)"{v) = zg{y) cj){y) for -l<y<l, 

In particular, one considers in this context the discrete string consisting of 
point masses connected by weightless thread, so that g is not a function but 
a linear combination of Dirac delta distributions. Then the solution to the 
inverse spectral problem can be expressed in terms of orthogonal polynomials 
and Stieltjes continued fractions [1, 2, 3, 26]. The reason for the appearance 
of Dirac deltas here is that when u has the form (1.8), the first derivative Ux 
has a jump of size —2mk at each point x — Xk, and this gives deltas in u^x 
when derivatives are taken in the sense of distributions. In each interval between 
these points, u is a linear combination of e^ and e~^, so Uxx = u there; thus 
m — u ~ Uxx — 2 X]fe=i "*fe ^Xk is ^ purely discrete distribution (or a discrete 
measure if one prefers). The measure g{y) in (1.9) is related to the measure m[x) 
through a so-called Liouville transformation, and g will be discrete when m is 
discrete. 

In the case of the Degasperis-Procesi and Novikov equations (and also for the 
Geng-Xue equation, as we shall see), the corresponding role is instead played by 
variants of a third-order nonselfadjoint spectral problem called the cubic string 
[22, 23, 20, 24, 17, 5]; in its basic form it reads 

-(j)"'{y)^ zg{y)(j){y) for -1 < y < 1, 
0(-l) = 0'(-l)=O, 0(1) = 0. 

The study of the discrete cubic string has prompted the development of a theory 
of Cauchy biorthogonal polynomials by Bertola, Gekhtman and Szmigielski 
[6, 5, 4, 7]; see Appendix A. In previous applications to peakon equations, the 
two measures a and (3 in the general setup of this theory have coincided (a = /3), 
but in this paper we will actually see two different spectral measures a and /3 
entering the picture in a very natural way. 

Like the above-mentioned PDEs, the Geng-Xue equation also admits peakon 
solutions, but now with two components. 



N 



{x,t) = y 'mk{t)e~ 


'\x-Xk(t)\ 


N 


\x-Xk{t)\ 

1 



(1.11) 



k=l 



where, for each k, at most one of rrik and n^ is nonzero (i.e., mkn^ = for all k). 
In this case, m and n will be discrete measures with disjoint support: 



N N 

m — u ~ u^_ 

k=l k=l 



xx = 2^mkSxk, n^v~Vxx^2^nkSx^. (1-12) 



This ansatz satisfies the PDE (1.2) if and only if the functions Xk{t), mk{t) and 
Ukit) satisfy the following system of ODEs: 

ik = u{xk)v{xk), 

rrik = mk{u{xk)vx{xk) - 2ux{xk)v{xk)) , (1-13) 

fi.k = nk{u.^{xk)v{xk) -2u{xk)v^{xk)), 
for k — 1,2, ... ,N. (Here we use the shorthand notation 

N 



and 



=1 



N 

Ux{xk) =^mi e-l^'=-^'l sgn(a;fe - x,). 



If u{x) = ^ TTii e~'^~^'', then the derivative u^ is undefined at the points Xk 
where rrik ^ 0, but here sgnO = by definition, so Ux{xk) really denotes the 
average of the one-sided (left and right) derivatives at those points. Note that 
the conditions rnk — and nik ^ both are preserved by the ODEs. Similar 
remarks apply to u, of course.) 

Knowing the solution of the inverse spectral problem for (1.3a)+(1.4a) in this 
discrete case makes it possible to explicitly determine the solutions to the peakon 
ODEs (1.13). Details about these peakon solutions and their dynamics will be 
published in a separate paper; here we will focus on the approximation-theoretical 
aspects of the inverse spectral problem. (But see Remark 4.12.) 

We will only deal with the special case where the discrete measures are 
interlacing, meaning that there are N = 2K sites 

Xi < X2 < ■ ■ ■ < X2K, 

with the measure m supported on the odd- numbered sites X2a-i, and the mea- 
sure n supported on the even- numbered sites X2a] see Figure 1 and Remark 3.1. 
The general formulas for recovering the positions Xk and the weights m2a-i and 
n2a are given in Corollary 4.5; they are written out more explicitly for illustration 
in Example 4.10 (the case K = 2) and Example 4.11 (the case K = 3). The case 
_?sr = 1 is somewhat degenerate, and is treated separately in Section 4.3. 
Appendix C contains an index of the notation used in this article. 

2 Forward spectral problem 

2.1 Transformation to a finite interval 

Let us start by giving a precise definition of the spectral problem to be studied. 
The time dependence in the two Lax pairs for the Geng-Xue equation will be 
of no interest to us in this paper, so we consider t as fixed and omit it in the 
notation. The equations which govern the x dependence in the two Lax pairs 
are (1.3a) and (1.4a), respectively. Consider the first of these: 



for xEK, (1.3a) 





where m{x) and n{x) are given. Our main interest lies in the discrete case, when 
m and n are actually not functions but discrete measures as in (1.12), but we 
will not specialize to that case until Section 3. 

There is a useful change of variables, similar to the one used for Novikov's 
equation [17], which produces a slightly simpler differential equation on a finite 
interval: 

y = tanha;, 
(t>i{y) = ipi (x) cosh X — ip^ (x) sinh x, 

(jjsiy) = z^Tp3{x)/coshx, (2.1) 

g{y) = ni{x) cosh x, 
h{y) — n{x) cosh'^ x, 

Under this transformation (with z 7^ 0), equation (1.3a) is equivalent to 

g{v) 02 , for -1< y < 1. (2.2a) 

/ V03/ 

(Notice that the 1 in the upper right corner of the matrix has been removed 
by the transformation. When h = g, equation (2.2a) reduces to the dual cubic 
string studied in [17].) In order to define a spectrum we impose the following 
boundary conditions on the differential equation (2.2a): 

M-^) = M-^)=0. 03(1) =0. (2.2b) 

By the eigenvalues of the problem (2.2) we then of course mean those values 
of A for which (2.2a) has nontrivial solutions satisfying (2.2b). 

The same transformation (2.1) applied to the twin Lax equation (1.4a) leads 
to the same equation except that g and h are interchanged. The spectrum of this 
twin equation will in general be different. To be explicit, the second spectrum is 
defined by the differential equation 





9{y) 








Hy) 


-A 







for -1 < y < 1, (2.3a) 

/ V^a/ 

again with boundary conditions 

02(-l)-03(-l)-O, 03(1) =0. (2.3b) 

Remark 2.1. Via the transformation (2.1), every concept pertaining to the 
original Lax equations (1.3a) and (1.4a) will have a counterpart in terms of the 
transformed equations (2.2a) and (2.3a), and vice versa. In the main text, we 
will work with (2.2a) and (2.3a) on the finite interval. However, a few things are 
more conveniently dealt with directly in terms of the original equations (1.3a) 
and (1.4a) on the real line; these are treated in Appendix B. More specifically, 
we prove there that the spectra defined above are real and simple, and we also 
obtain expressions for certain quantities that will be constants of motion for the 
Geng-Xue peakon dynamics. 



Remark 2.2. Transforming the boundary conditions (2.2b) back to the real 
line via (2.1) yields 

lim 7/12(2;) = lim e^'ipsix) = 0, lim e'^'tf^six) = 0. (2.4) 

2;— >- — CO a;— ^ — C30 cc— ^+00 

Each eigenvalue A 7^ of (2.2) corresponds to a pair of eigenvalues z = ±\/— A 
of (1.3a)+(2.4). As an exceptional case, A = is an eigenvalue of (2.2), but 
z = is not an eigenvalue of (1.3a)+(2.4); this is an artifact caused by the 
transformation (2.1) being singular for z — 0. When talking about eigenvalues 
below, we will refer to A rather than z. 

In Section 2.3 below we will also encounter the condition 0i(— 1) = 1; this 
translates into 

lim e-%M^; z) - M^; z)) ^ 2. (2.5) 

2:— > — C30 

2.2 Transition matrices 

Let 

h{y) 
^(y;A)=| g{y) \ , A{y; X) ^ \ %) | (2.6) 

^A 

denote the coefficient matrices appearing in the spectral problems (2.2) and 
(2.3), respectively. To improve readability, we will often omit the dependence 
on y in the notation, and write the differential equations simply as 

_=^(A)ci>, ^ = Aim (2.7) 

respectively, where $ = {(pi, (j)2, ^s)"^- Plenty of information about this pair of 
equations can be deduced from the following modest observation: 

Lemma 2.3. The matrices A and A satisfy 

^0 1\ 
AiX) = -JA{~Xyj, where J=|0 -1 = J^ = J"\ (2.8) 

0/ 

Proof. A one-line calculation. D 

Definition 2.4 (Involution a). Let a denote the following operation on the (loop) 
group of invertible complex 3x3 matrices A" (A) depending on the parameter 
Ae C: 

A(A)'" = JA(-A)-^ J. (2.9) 

(We use the customary abbreviation X^^ = (A-^)^^ = (A~^)^.) 

Remark 2.5. It is easily verified that a is a group homomorphism and an 
involution: 

(A(A)r(A))^ = x{xrY{xr, {x{xry = a(a). 



Definition 2.6 (Fundamental matrices and transition matrices). Let U{y; A) 
be the fundamental matrix of (2.2a) and U{y; A) its counterpart for (2.3a); i.e., 
they are the unique solutions of the matrix ODEs 

dU 

— ^A{y;X)U, C/(-l;A)=J, (2.10) 

and _ 

^=^(y;A)t/, C7(-1;A) = /, (2.11) 

oy 

respectively, where / is the 3x3 identity matrix. The fundamental matrices 
evaluated at the right endpoint y ~ I will be called the transition matrices and 
denoted by _ _ 

Six) = L/(l; A), S{X) = C/(l; A). (2.12) 

Remark 2.7. The fundamental matrix contains the solution of any initial value 
problem: $(?/) = U{y; A)$(— 1) is the unique solution to the ODE d^/dy — 
A{X)^ satisfying given initial data $(—1) at the left endpoint y — —1. In 
particular, the value of the solution at the right endpoint j/=lis$(l) = 

siXM-i). 

Theorem 2.8. For all y G [-1, 1], 

det [/(y; A) = det C/(y; A) = 1 (2.13) 

and 

U{y-X) = U{y-XY. (2.14) 

In particular, det <S'(A) — det 5(A) = 1, and S{X) = S{XY . 

Proof. Equation (2.13) follows from Liouville's formula, since A is trace-free: 

j-y 
det U{y; X) = (det C/(-l; A)) exp / tr A{£,; X) di = (det /) exp = 1, 

and similarly for U. To prove (2.14), note first that 

dy dy 

= -C/(A)-i A{X) t/(A) t/(A)-i = -C/(A)-i AiX), 

which implies that 



-U U[-X) " , 

/ac/(-A)-i^^ 



dy dy 



= J{-U{-X)-^A{-X)fj 

= -jA{-xfu{-xy^j. 



With the help of Lemma 2.3 this becomes 

^u{xy = A{x)Ju{-x)'^j 
^A{x)u{xy. 

Since U{Xy = I = U{X) when y = -1, we see that U{Xy and U{X) satisfy 
the same ODE and the same initial condition; hence they are equal for all y by 
uniqueness. D 

Corollary 2.9. The transition matrices 5(A) and S{X) satisfy JS{X)^ JS{—X) = 

I and S{—X) = J(adj S{X)) J, where adj denotes the adjugate (cofactor) matrix. 
In detail, this means that 





(2.15) 



and 



S11S22 — S12S21 S11S23 — S21S13 S12S23 ~ <S'i3S'22 

'S'( — A) = I 511532 — 5i253i 5ii533 — 5i353i 5i2533 — 5i3532 I . (2.16) 

^52i532 — 522531 52i533 — 523531 522533 — 523532y 

(The subscripts ±A indicate the point where the matrix entries are evaluated.) 

2.3 Weyl functions 

Consider next the boundary conditions 02(~1) = </'3(— 1) = = (^3(1) in 
the two spectral problems (2.2) and (2.3). Fix some value of A G C, and let 
$ — {(j)i,(j)2,4>3)'^ be a solution of d^/dy — A{X)<^ satisfying the boundary 
conditions at the left endpoint: (?!)2(— 1) = ^3(^1) — 0. For normalization, 
we can take (/ii(— 1) — 1; then the solution $ is unique, and its value at the 
right endpoint is given by the first column of the transition matrix: $(1) = 
5(A)$(-1) = 5(A)(1,0,0)'^ = (5ii(A), 521(A), 53i(A))^. This shows that the 
boundary condition at the right endpoint, (jiz^i) = 0, is equivalent to 531(A) — 0. 
In other words: A is an eigenvalue of the first spectral problem (2.2) if and only 
if 531(A) = 0. 

We define the following two Weyl functions for the first spectral problem 
using the entries from the first column of 5(A): 

The entries of 5(A) depend analytically on the parameter A, so the Weyl functions 
will be meromorphic, with poles (or possibly removable singularities) at the 
eigenvalues. The signs in (2.17) (and also in (2.18), (2.24), (2.26) below) are 
chosen so that the residues at these poles will be positive when g and h are 
positive; see in particular Theorem 3.10. 



Similarly, A is an eigenvalue of the twin spectral problem (2.3) if and only if 
'S'3i(A) = 0, and we define corresponding Weyl functions 

W{X) = ~^^, Z{X)^~^^. (2.18) 

^3l(A) 53l(A) 

Theorem 2.10. The Weyl functions satisfy the relation 

Z{X) + WiX)Wi-X) + Z{-X)=0. (2.19) 

Proof. The (3, 1) entry in the matrix equality (2.15) is 

^3i(A)5ii(-A) - S2i{X)S2i{~X) + Sn{X)S3i{~X) = 0. 
Division by — S'3i(A)S'3i(— A) gives the desired result. D 

2.4 Adjoint spectral problems 

Let us define a bilinear form on vector- valued functions 




(2.20) 



with0i,(/)2>3 €^^(-1,1): 

{^,n)^J ^yfjn{y)dy 

{<j)i{y)uj3{y) ~ (j)2{y)uj2{y) + cl)3{y)^i{y)) dy. 
Lemma 2.3 implies that 

('a(a)$) jn = -<i>'^j(A{-x)n), 

which, together with an integration by parts, leads to 

{{i~MX)) 'J''") - ["^"^jK^-i - {'^^ (I; -^(-^)) ">• (2-21) 

Now, if $ satisfies the boundary conditions (/)2(— 1) = (^^{—l) = = 03(1), then 
what remains of the boundary term [<i>-^J5l] is 

01(1)W3(1) - 02(1)^2(1) - (/)l(-l)cc-3(-l), 

and this can be killed by imposing the conditions W3(— 1) = = 102(1) = ^3(1). 
Consequently, when acting on differentiable L^ functions with such boundary 
conditions, respectively, the operators -j — ^(A) and — ^ + A{—X) are adjoint 
to each other with respect to the bilinear form (■,■): 



This calculation motivates the following definition. 



Definition 2.11. The adjoint problem to the spectral problem (2.2) is 

r)fl ~ 

— = A{-X)n, (2.22a) 
dy 

W3(-1) = = W2(1) = W3(1). (2.22b) 

Proposition 2.12. Let ^{y) be the unique solution of (2.22a) which satisfies 
the boundary conditions at the right endpoint, W2(l) — ^3(1) — 0, together with 
wi(l) — 1 (for normalization) . Then, at the left endpoint y — —1, we have 

/1\ /1\ /^33(A)\ 

f}(-l) = S{-X)-' U JS{Xfj p U -532(A) . (2.23) 

\0j \0j V ^3i(A) J 

Proof. Since (2.22a) agrees with the twin ODE (2.3a) except for the sign of A, 
the twin transition matrix with A negated, S{—X), will relate boundary values 
of (2.22a) at y = -1 to boundary values a.t y = +1: (1,0,0)^ = ^{1) = 
S{-X)n{-1). The rest follows from Corollary 2.9. D 

Corollary 2.13. The adjoint problem (2.22) has the same spectrum as (2.2). 

Proof. By (2.23), the remaining boundary condition W3(— 1) = for (2.22) is 
equivalent to 531(A) = 0, which, as we saw in the previous section, is also the 
condition for A to be an eigenvalue of (2.2). D 

We define Weyl functions for the adjoint problem as follows, using the entries 
from the third row of S{X) appearing in (2.23): 

'-'3l(,Aj J3U'\j 

To complete the picture, we note that there is of course also an adjoint 
problem for the twin spectral problem (2.3), namely 

— = A{~X)fl, (2.25a) 
dy 

W3(-l) = = 0.2(1) = ^3(1). (2.25b) 

A similar calculation as above shows that the eigenvalues are given by the zeros 
of 531 (A), and hence they are the same as for (2.3). We define the twin adjoint 
Weyl functions as 

W*{X)^-^^, Z*{X) = -&^. (2.26) 

531(A)' 531(A) 

Theorem 2.14. The adjoint Weyl functions satisfy the relation 

Z*{X) + W*{X)W*{-X) + Z*{-X) = 0. (2.27) 

Proof. Since a matrix commutes with its inverse, we can equally well multiply 
the factors in (2.15) in the opposite order: 5(— A) • JS{X)^J = I. Division of 
the (3, 1) entry in this identity by — 53i(A)53i(— A) gives the result. D 



10 



3 The discrete case 

We now turn to the discrete case (1.12), when m(x) and n{x) are discrete 
measures (linear combinations of Dirac deltas) with disjoint supports. More 
specifically, we will study the interlacing discrete case where there are N = 2K 
sites numbered in ascending order, 

xi < a;2 < • •• < X2K, 

with the measure m supported on the odd-numbered sites X2a-i, Siiid the 
measure n supported on the even- numbered sites X2a- That is, we take TO2 = 
7714 = • • ■ = and ni = 713 = • • • = 0, so that 



N K 

m 

k=l a=l 

N K 

n 

k=l a=l 



2 2J "7i; (5a;;, ^ 2 zJ ™2a-l Sx2a-i : 

(3.1) 
2^^nkSx^ = 2 ^ n2a ^x^^ ■ 



We will also assume that the nonzero nik and Uk are positive; this will be needed 
in order to prove that the eigenvalues A — —z^ are positive. The setup is 
illustrated in Figure 1. 

Given such a configuration, consisting of the AK numbers {xk^ rn2a-i, n2a}, 
we are going to define a set of AK spectral variables, consisting of 2K — 1 eigen- 
values Ai, . . . , Xk and fii, . . . , hk-i, together with 2K + 1 residues oi, . . . , qk, 
bi,. . . ,6/f_i, 6oo and 6^. In Section 4 we will show that this correspondence 
is a bijection onto the set of spectral variables with simple positive ordered 
eigenvalues and positive residues (Theorem 4.8), and give explicit formulas for 
the inverse map (Corollary 4.5). 

Remark 3.1. Non-interlacing cases can be reduced to the interlacing case by 
introducing auxiliary weights at additional sites so that the problem becomes 
interlacing, and then letting these weights tend to zero in the solution of the 
interlacing inverse problem; the details will be published in another paper. 

Remark 3.2. The case K — I is somewhat degenerate, and also rather trivial. 
It is dealt with separately in Section 4.3. In what follows, we will (mostly without 
comment) assume that K > 2 whenever that is needed in order to make sense of 
the formulas. 

Under the transformation (2.1), the discrete measures m and n on R are 
mapped into discrete measures g and h, respectively, supported at the points 

Uk = tanh.-Efc (3.2) 

in the finite interval (—1,1). The formulas g{y) = 7n-(a;) cosh^ a; and h{y) — 
n{x) cosh^ X from (2.1) should be interpreted using the relation 6{x — Xkjdx = 
5{v - yk)dy, leading to Sx^{x) = Sy^{y)^{xk) = Sy^{y) / cosh^ Xk- Since we 
will be working a lot with these measures, it will be convenient to change the 
numbering a little, and call the weights gi,g2T ■ ■ ,9k and hi,h2, ■ . ■ , hj^ rather 
than (7i, (73, . . . , g2K-i and /12, /14, • • ■ , h2K', see Figure 2. With this numbering, 
we get 

K K 

9 = ^9aSy2^_^, h = ^haSy^^, (3.3) 

a=l a—1 
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2mi 



2n2K 



2m3 



2n4 



2n2 



2m2K-i 



Xi X2 



X3 



X4 



X2K-1 



X2K X 



Figure 1. Notation for the measures m and n on the real line R in the interlacing 
discrete case (3.1). 



91 



hi 



52 



QK 



flK 



_i yi 2/2 

< — >< — >< — 

lo h I2 



Vi 



Vi 



h 



y2K-i 



h 



y2K 



h 



+1 



Figure 2. Notation for the measures g and h on the finite interval ( — 1, 1) in the 
interlacing discrete case (3.3). 



where 



5a = 2m2a-i coshxaa-i, 



2n2aCOsha;2a- 



(3.4) 



3.1 The first spectral problem 

The ODE (2.2a), dy^ = A{y; X)^, is easily solved explicitly in the discrete case. 
Since g and h are zero between the points y^, the ODE reduces to 




in those intervals; that is, (f)i and (f)2 are constant in each interval yk < y < yu+i 
(for < A: < 2K, where we let 2/0 = ^1 and j/2A'+i = +1), while (/)3 is piecewise 
a polynomial in y of degree one. The total change in the value of (j)^ over the 
interval is given by the product of the length of the interval, denoted 



h 



yk+i -yk, 



(3.5) 



12 



and the slope of the graph of (f)^; this slope is —A times the constant value of 0i 
in the interval. In other words: 



where the propagation matrix Lk is defined by 




(3.6) 



(3.7) 



At the points y^, the ODE forces the derivative 9j,$ to contain a Dirac delta, and 
this imposes jmxip conditions on $. These jump conditions will be of different 
type depending on whether k is even or odd, since that affects whether the Dirac 
delta is encountered in entry (1,2) or (2,3) in the coefficient matrix 



h{y) 

9{y) 




A{y;\) = 
More precisely, when fc = 2a is even, we get a jump condition of the form 

Hyt)-Hvu) 



fQ ha 0^ 

I $(2/fe). 

^0 0; 



This implies that (^2 and cf)^ don't jump at even-numbered yk, and the continuity 

of 4)2 in particular implies that the jump in 0i has a well-defined value ha<j)2{y2a)- 

When fc = 2a — 1 is odd, the condition is 

/O 0' 

Hyt)-Hvk)^ P 9a\Hyk). 
Vo 0, 

Thus, (pi and 03 are continuous at odd-numbered j/fc, and the jump in (/)2 has a 
well-defined value ga4>3{y2a-i)- 

This step-by-step construction of $(y) is illustrated in Figure 3 when $(—1) = 
(1,0, 0)"^; as we have already seen, this particular case is of interest in connection 
with the spectral problem (2.2) where we have boundary conditions 02(^1) = 

03(-l) =0 = 03(1). 



With the notation 



1 


X 


^xy 





1 


y 








1 



(3.8) 



the jump conditions take the form 



^(yl) 



hr, 



Hy2a)^ ^(yL-i) 



*(y2a-l)- 



(For the purposes of this paper, the top right entry of [^1 might as well have 
been set equal to zero; we have defined it as ^^y only to make [^] agree with 



13 



(?ii(y;A) (solid) 



t 



yi 



y2 



ys 



yi 



ys 



ye 



4>2 (y; A) (dashed) 



'■■•■.. ...•■■■' 03(y;A) (dotted) 

Figure 3. Structure of the solution to the initial value problem dy^ — A{y; A)$ 
with $(— 1; A) = (1,0,0)"^, in the discrete interlacing case. The components (f>i 
and (ji2 are piecewise constant, while cjf>3 is continuous and piecewise linear, with 
slope equal to —A times the value of cpi. At the odd-numbered sites j/2a-i, the 
value of <j!>2 jumps by (?a<?f>3(2/2a-i)- At the even-numbered sites y2a, the value of 
01 jumps by ha4>2{y2a)- The parameter A is an eigenvalue of the spectral problem 
(2.2) iff it is a zero of 03(1; A), which is a polynomial in A of degree K + 1, with 
constant term zero. This picture illustrates a case where A and the weights ga 
and ha are all positive. 



a jump matrix appearing in our earlier work [23, 17].) We can thus write the 
transition matrix ^(A) as a product of 1 -|- 4:K factors, 



5(A) = L2kW 







L2K-l{^) 



■ ■ 

.9K, 



L2K-2{^) 



Li(A) 



LoW- (3.9) 



For later use in connection with the inverse problem, we also consider the partial 
products Tj{X) containing the leftmost 1 -|- 4j factors (for j — 0, . . . ,K); put 
differently, Tk-j{X) is obtained by omitting all factors after L2j{X) in the product 
for 5(A): 



TK-j{X)^L2Ki\)' 







^2j + l 



(A) 





.5j+i 



L2,(A). 



(3.10) 



Thus 5(A) = Tk{X), and Tk-;j{X) depends on {gj+i, . . . .qk), (/ij+i, . . . , /lif), 
(/2j, . . • ,?2i<-)- 

Proposition 3.3. The entries ofTj(X) are polynomials in X, with degrees as 
follows: 

( i j-i j-i^ 

degr,(A)- J j-1 j-l| (j>l). (3.11) 

\i + 1 3 i 



The constant term in each entry is given by 



Tk-M 







■ 0' 

9K, 



T-j + l 





9]+i 



a>j O'^^yj 

1 J29a 



Vo 



a>j 



(3.12) 
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For those entries whose constant term is zero, the coefficient of X^ is given by 

/ * * *\ 

2a-2 
dTK-j 



dX 



-(0) 



V 



* 


* 


2a-2 




> ! > ; 9alk 


* 


a>j k=2j 




2K 


2K 


->:i. 


->:>>jk 


k=2j 


a>j k—2a 



(3.13) 



/ 



The highest coefficient in the (3, 1) entry is given by 

(T,)3l(A) = (-AV + W n h„A I n 9aha 

\m—K—j j \a— K+1— J 



(3.14) 



Proof. Equation (3.12) follows at once from setting A = in (3.10). Next, group 
the factors in fours (except for the lone first factor L2k{X)) so that (3.10) takes 
the form Tk-j = L2K tK tx-i ■ ■ ■ tj+i, where 



ta{X) 



L2a-l{X) 



i2a-2(A)=(0 1 3„ -A b.-2g„ (1,0,0). 

^0 1 / V'2a-l+i2a-2/ 



The degree count (3.11) follows easily by considering the highest power of A arising 
from multiplying these factors, and (3.14) also falls out of this. Differentiating 
Tr-_j+i(A) = Tk-j{X) tj{X) and letting A = gives 



dT, 



K-j + l 



(0) 



"^ ""-' (0) 1 gf - Tk^,{0) h,-29, 0. (3.15) 

VO 1 / \l2,-l+l2j-2 0^ 



dX ^ ' dX 
With the help of (3.12) one sees that the (3, 1) entry of this equality reads 

{T'K^,+i)3m - {T'K^.hm - ihj-i + 12,-2), (3.16) 

the (3,2) entry is 

(r^„,+i)32(0) = (t;,_^.)32(o) + h,iTK-,y,M, (3.17) 

and the (2, 1) entry is 

{T'k-,+i)2i{0) = (r;._^.)2i(0) - ;2,-25, - (T.9a) (%-i + ^2,-2). (3.18) 

Solving these recurrences, with the inital conditions coming from T'q(O) — L'^j^i^) 
(i.e., —I2K in the (3, 1) position, zero elsewhere), gives equation (3.13). D 

We also state the result for the important special case S{X) = Tk{X). (In 
the (3, 1) entry, J2k=o 'fe = 2 is the length of the whole interval [—1, 1].) 

Corollary 3.4. The entries of S{X) are polynomials in X, with 



K K-l K-V 
deg5(A)= ( K K-l K-l 
,K + l K K 



(3.19) 
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5(0) 



a a>b 

1 E^a 

a 

\o 1 / 



(3.20) 



5'(0) 



/ * * 

/^ 2a-2 

~ Y. J2 9Jk * 

a=l fc=0 

2^" K 2_ff 

\ fc=0 a=l fc=2a 

53l(A) = (-A)^'^'+l ( n hm] ( n .ga/ia ) 

\m=0 / \a=l / 



*\ 



(3.21) 



(3.22) 



Remark 3.5. The ideas that we use go back to Stieltjes's memoir on continued 
fractions [29] and its relation to an inhomogeneous string problem, especially its 
inverse problem, discovered by Krein in the 1950s. A comprehensive account 
of the inverse string problem can be found in [11], especially Section 5.9. The 
connection to Stieltjes continued fractions is explained in [13, Supplement II] 
and in [1]. Briefly stated, if (j){y] A) satisfies the string equation 



(t>yy = >'9{y)(t>i 



-l<y< 1, 



-1;A)=0, 



with a discrete mass distribution g{y) — Yl^=i9j^yj^ then the Weyl function 
W{\) = I(~i.\\ admits the continued fraction expansion 



W{z) 



ln + 



-zg-n 



ln~l + 



+ 



-zg2 



h + 



-zgi 



lo 



Q2,(A) 



(where Ij — Vj+i — Hj), whose convergents (Pade approximants) T2j{\) = ^ 
satisfy 

/ n—l 

Py{X) = {-iygn( n ^kgk)x' + 
Q2,{X) = {-iy[ n hgk)x' + 

^k—n—j-\-l 
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(t>i{y;X) (solid) 



-l'\ 



yi y2 y-i va y^ ye 



4>'i{y\ A) (dashed) 



03 (yi -^) (dotted) 



Figure 4. Structure of the solution to the twin problem 9y<E> = .4(y; A)$ with 
$(— 1; A) = (1,0,0)'^, in the discrete interlacing case. The differences compared 
to Figure 3 are the following: At the odd- numbered sites y2a-i, the value of (j>i 
(not (j)2) jumps by ga4'2{y2a-i) ■ At the even-numbered sites j/2a, the value of (j>2 
(not 0i) jumps by ha<j>3{y2a)- The parameter A is an eigenvalue of the twin spectral 
problem (2.3) iff it is a zero of (^3(1; A), which is a polynomial in A of degree K 
(not K + 1), with constant term zero. Note that the first mass gi has no influence 
here. (Indeed, since 4>2{y\', A) = 0, there is no jump in (j>\ at y = j/i, regardless of 
the value of (ji.) 



3.2 The second spectral problem 

For the twin ODE (2.3a), dy^ — A{y\ A)$, where the measures g and h are 
swapped, the construction is similar. The only difference is that the weights 
Qa at the odd-numbered sites will occur in the type of jump condition that we 
previously had for the weights ha at the even- numbered sites (and vice versa). 
Thus, the transition matrix is in this case 



S{\) = L2k{X) 



■ 



L2K-l{^) 



9K 





L2K-2{^) 



Li(A) 



Lo(A). (3.23) 



This solution is illustrated in Figure 4 for the initial condition $(—1) = (1,0,0)"^. 
It is clear that it behaves a bit differently, since the first weight gi has no 
influence on this solution $, and therefore not on the second spectrum either. 
(The first column in [^^]Lo(A) does not depend on gi.) 

Let Tj{X) be the partial product containing the first 1 + 4j factors in the 
product for S{X); in other words. 



Tk-,{X)^L2kW^ 







tj+i 



^2j + l 



(A) 



93 + 1 





L2,W- 



(3.24) 



Proposition 3.6. The entries ofTj(X) are polynomials in X, satisfying 

1 9K 

ri(A)=( -XhK{l2K-l + l2K-2) l-^hKgKhK-i hx \ , (3.25) 

-X{l2K + hK~i + hK-2) —XgxihK + I2K-1) 



degT,(A) 



i-1 


j-i J -2 


J 


J j-i 


J 


J j - 1 



(j > 2), 



(3.26) 
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Tk-A^) 



' ■ 



9k 








9]+i 




(^ ^9a ^ 9ahb\ 

a>j a>b>j 

1 ^/la 

Vo 



a>j 



(3.27) 



/ 



dX 



(0) 



* 

2a- 1 







^\ 



(3.28) 



a>j k=2j 
2K 

k=2j 

For < j < K — 2, the highest coefficients in the (2, 1) and (3, 1) entries are 
given by 



V 



2K 

"E E 9alk * 
a>j fc— 2a— 1 



(Tk_,)2i(A) = (-A)^-^ I fl hm-l I (^2, + l + ^2^) hK ( J] 3«+l^« 

, m=j+2 



i-xy^-^hK + hK~i) n ^2m-i I (%+i + hj) ( n 5«+l^« 

\ m=j+2 



(3.29) 



(3.30) 



where Jlm=i<' ~ ^- Moreover, Tj and Tj are related by the involution a (see 
Definition 2.4): 

f,(A) = T,{Xr. (3.31) 

Proof. The degree count and the coefficients are obtained Hke in the proof of 
Proposition 3.3, although the details are a bit more involved in this case. (Group 
the factors in Tj as follows: L2i<"(A) times a pair of factors, times a number a 
quadruples of the same form as ta{X) in the proof of Proposition 3.3 but with 
ha and ga replaced by ga+i and ha respectively, times a final pair at the end.) 
The CT-relation (3.31) can be seen as yet another manifestation of Theorem 2.8, 
and (since cr is a group homomorphism) it also follows directly from the easily 
verified formulas LkiX)" = Lk{X) and g] " = [^] . □ 

We record the results in particular for the case S{X) = Tfj-(A): 

Corollary 3.7. The entries of S{X) are polynomials in A, satisfying 

^K-1 K-1 K-2'^ 
degS'(A)= \ K K i^ - 1 I , (3.32) 



K 



K K~ 1 
K K- I, 



/I E^" E^^M 



a>b 



5(0) = 





Vo 



a 

1 / 



(3.33) 
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5'(0) = 



/ * 


* * 


K 2a- 1 




->;>;w. 


* * 


a=l fe=0 






K 2K 


-2 


~y, y, 9ah * 

a=l k=2a-l 



(3.34) 



('T/ie interpretation of (3.32) when K = 1 is that the (1,3) entry is the zero 
polynomial.) The leading terms 0/^21 (A) and 6*31 (A) are given by 



(3.35) 



521(A) = (-A)^' [ n hm-l ) {lo + h) hK ( n 9a+lha\ + ' ' - 
\m=2 / \a=l / 

§3l(A) = (-A)^(/2K + I2K-1) ( n ^2m-l ) Go + '1) f n 5a+l/la ) + ' ' ' , 

\m=2 / \a=l / 

(3.36) 

with the exception of the case K — 1 where we simply have 5*31 (A) = — 2A. (The 
empty product nm=2 ^2m-i is omitted from S'31 in the case K = 2, and from S21 
in the case K ^ I.) 

3.3 Weyl functions and spectral measures 

Since the entries of 5(A) are polynomials, the Weyl functions W = —S21/SS1 
and Z = —S11/S31 are rational functions in the discrete case. They have 
poles at the eigenvalues of the spectral problem (2.2). Likewise, the twin Weyl 
functions W = —S21/S31, Z = —S11/S31 are rational functions, with poles at 
the eigenvalues of the twin spectral problem (2.3). 

Theorem 3.8. // all g^ and hk are positive, then both spectra are nonnegative 
and simple. The eigenvalues of (2.2) and (2.3) will be denoted by 



= Aq < Ai < • • • < A 



K 



= Ha < Hi < ■ ■ ■ < fiK-i 



(zeros of S31), 
(zeros of S31). 



(3.37) 
(3.38) 



Proof. This is proved in the appendix; see Theorem B.l. (It is clear that if 
the zeros of the polynomials S'31 (A) and S'31 (A) are real, then they can't be 
negative, since the coefficients in the polynomials have alternating signs and 
all terms therefore have the same sign if A < 0. However, it's far from obvious 
that the zeros are real, much less simple. These facts follow from properties of 
oscillatory matrices, belonging to the beautiful theory of oscillatory kernels due 
to Gantmacher and Krein; see [13, Ch. II].) D 

Remark 3.9. Taking Propositions 3.4 and 3.7 into account, we can thus write 



K 



531(A) =-2Ani-^ 



i=l 



S3i(A) 



-2A J] ( 1 



A 



(3.39) 
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Theorem 3.10. // all g^ and hk are positive, then the Weyl functions have 
partial fraction decompositions 

K 



. A - A, 

■l—l 



^-^ b 



W{\) ^-hoo+Y. Y^' (^-"^^^^ 



K 



H 



^(A) = ^ + ElAr^ (3.40c) 



2A ^^ A - A, 

1 ^^ d, 



^(^) = ^ + Eir^' (3.40d) 

2A ■'^-^ A — u, 

i=i -^ 

where Oj, 6j, 6oo, Q, dj are positive, and where W and W determine Z and Z 
through the relations 

j — l •' i—l •' 

Proof The form of the decompositions foUows from Propositions 3.4 and 3.7 
(polynomial degrees), together with Theorem 3.8 (all poles are simple). In 
W ~ —S21/S31 the factor A cancels, so there is no residue at A = 0, and similarly 
for W — —S21/S31 (which however is different from W in that the degree of the 
numerator equals the degree of the denominator; hence the constant term — 6oo)- 
The residue of Z =^-811/831 at A = is -S'ii(0)/5'^i(0) = 1/2 by Corollary 3.4, 
and similarly for Z(X). 

From the expressions (3.35) and (3.36) for the highest coefficient of 5*21 
and 5*31 we obtain (for K > 2) 

&00 = - lim WiX)= lim |ii^ = J^EhE^ , (3.42) 

which shows that 600 > 0. (In the exceptional case K — 1 we have instead 
-W{X) = lhi{l2 + li) = boo>0.) 

The proof that a^ and bj are positive will be given at the end of Section 3.5. 
It will then follow from (3.41) that Ci and dj are positive as well. 

To prove (3.41), recall the relation Z{X) + W{X)W{-X) + Z{-X) = from 
Theorem 2.10. Taking the residue at A = Ai on both sides yields 

c, + a,W{~X,) + = 0. 

Taking instead the residue at A = fij in Z{-X) + W{-X)W{X) + Z{X) = 0, we 
obtain 

+ W{-^lJ)bj + dj =0. 

n 
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Definition 3.11 (Spectral measures). Let a and /3 denote the discrete measures 

K K-\ 

a = ^ a,;(5A, , /^ = X! ^i^i^i ' (3.43) 

where a^ and hj are the residues in W{X) and Ty(A) from (3.40a) and (3.40b). 

We can write W and W in terms of these spectral measures a and /3, and 
likewise for Z and Z if we use (3.41): 

inA) = /f*^. («4.) 

(Note the appearance here of the Cauchy kernel l/(x + j/).) 

We have now completed the spectral characterization of the boundary value 
problems (2.2a) and (2.3a). The remainder of Section 3 is devoted to establishing 
some basic facts which will be needed for formulating and solving the inverse 
problem in Section 4. 

3.4 Rational approximations to the Weyl functions 

The Weyl functions W(X) and Z(X) are defined using entries of the transition 
matrix S(X). Next, we will see how entries of the matrices Tj{X) (partial products 
of S{X); see (3.10)) produce rational approximations to the Weyl functions. We 
have chosen here to work with the second column of Tj{X), since it seems to 
be the most convenient for the inverse problem, but this choice is by no means 
unique; many other similar approximation results could be derived. 

Theorem 3.12. Fix some j with 1 < j < K , write T{\) ~ Tj(X) for simplicity, 
and consider the polynomials 

Q(A) = -T32(A), P(A)=T22(A), i?(A) = ri2(A). (3.45) 

Then the following properties hold: 

degQ^j, degP = j-l, degi? = j-l, (3.46) 

g(0) = 0, P(0) = 1, (3.47) 

and, as X -^ oo, 

W{X)QiX) - PiX) = O (jY (3.48a) 

Z{X)Q{X) - R{X) ^ O (jY (3.48b) 

R{X) + P{X)W{-X) + Q{X)Z{-X) - O (^) . (3.48c) 

(For j = K , the right-hand side of (3.48c) can he replaced by zero.) 
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Proof. Equations (3.46) and (3.47) were already proved in Proposition 3.3. With 
the notation used in that proof, the first column of the transition matrix S{X) is 
given by 

%i(A)\ 

S21W = L2k{\) tK{\) ■ ■ ■ tK+l^j{\) tK-j{\) ■ ■ ■ h{X) ( ) , 

.^31 (A)/ ^ '■ ^"^■ 



=T{\) 



where ai, 02, 03 have degree at most K — j in A. Hence, 

•JSl 'J31 



Qi(A) 
a2(A) 
Q3(A) 



J32(T2i, T22, r23) ( °^ j - r22(T3i, T32, T33) (02! 













S31 




-fll 


T21 
T31 


T22 
T32 


+ 


as 


T22 
T32 


T23 
T33 



-ai(T-i)3i+a3(T-i) 
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'5'31 



'S'31 



where the last step uses that det T(A) = 1 (since each factor in T has determinant 
one). By (3.31), T-i(A) = jf(-A)'^J, where r(A) is shorthand for T,(A) 
(defined by (3.24)). In particular, (r-i)3i(A) = T3i(-A) and (T-i)ii(A) = 
r33(-A), so 



M/(A)Q(A) - P(A) 
By (3.19) and (3.26) we have 



-ai(A)r3i(-A)+a3(A)r33(^A) 
531(A) 



deg53i-i^+l, degT3i=j, degT33=j-l, 

which shows that WQ - P = [\(^-f'>+:>-(^+^)) = O {\-^) as A ^ 00. 
The proof that ZQ — R = O (A^^) is entirely similar. 
To prove (3.48c), we start from 

^21 (A) =r(A) 62(A) 
,^31 (A)/ \hW, 

where 61, 62, 63 have degree at most K — j. Using again T(A) = T{XY = 
JT{-\)-^J, we obtain 

/Sii(-A)\ 

-62(-A) = -(0,l,0)jr(A)^J S2i(-A) 

\S3i(-A)/ 

/ 2;u(A) \ 

= (0,1,0)T(A)^ -5,i(A) 

\Sii(-A)/ 

= (i?(A),P(A),-Q(A))f wix) \s3i{\). 

\-Z(-X)J 
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Since S31 has degree K by (3.32), we find that R{X)+P{X)W{-X)+Q{X)Z{-X) = 
-62(-A)/S'3i(A) = O (A^^-^)-^) = O (A-J). (When j ^ K we have 62(A) - 
0.) D 

Remark 3.13. Loosely speaking, the approximation conditions (3.48) say that 

PW ^ wfA^ ^(^) - Z(X) 
Q(A)-^(")' Q(A)-^(')' 

and moreover these approximate Weyl functions satisfy 

|(A) + §(A)M?(-A) + Z(-A)«0 
in place of the exact relation 

Z{X) + W^(A)W^(-A) + Z{-X) = 

from Theorem 2.10. We say that the triple (Q, P, R) provides a Type I Hermite- 
Pade approximation of the fmictions W and Z^and simultaneously a Type II 
Hermite-Pade approximation of the functions W and Z; see Section 5 in [6]. 

We will see in Section 4 that for given Weyl functions and a given order of 
approximation j, the properties in Theorem 3.12 are enough to determine the 
polynomials Q, P, R uniquely. This is the key to the inverse problem, together 
with the following simple proposition. We will need to consider Q, P, R for 
different values of j, and we will write Qj: Pjj Rj to indicate this. As a somewhat 
degenerate case not covered by Theorem 3.12 (the degree count (3.46) fails), we 
have 

Qo(A) = 0, Po(A) = l, i?o(A) = 0, (3.49) 

coming from the second column of To (A) = L2k{X). 

Proposition 3.14. // all Qj{X) and Rj{X) are known, then the weights hj and 
their positions y2j can be determined: 

hj=RK-j+m-RK-M. (3.50) 

(1 ~ y2j)h, - Q'k-,+M - Q'^_,(0), (3.51) 

/or j = 1, . . . , K. 

Proof. By definition, Qj = — (Tj)32 and Rj = (Tj) 12, and Proposition 3.3 says 
that Rk-M = Ea>, ha and Q'j,_,{0) = Ea>, Efefsa hJk, ioTO<j<K-l. 
The statement follows. (Note that J2k=2i ^k — ^ ~ V2j-) D 

In order to access the weights gj and their positions 1/2 j-i we will exploit the 
symmetry of the setup, via the adjoint problem; see Section 3.6. 

3.5 Adjoint Weyl functions 

Recall the adjoint Weyl functions defined by (2.24) and (2.26), 

W* ^-832/831, Z* = -833/831, W* = -832/831, Z* = -833/831, 
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which have the same denominators as the ordinary Weyl functions 

W^-S2i/S^i, Z^~Sii/S3i, W =-821/831, Z =-811/831, 

but different numerators. Since the transition matrices 8{X) and 8{X) both have 
the property that the (2,1) and (3,2) entries have the same degree, and the 
(1, 1) and (3,3) entries have the same degree (see Propositions 3.4 and 3.7), the 
adjoint Weyl functions will have partial fraction decompositions of exactly the 
same form as their non-starred counterparts (cf. Theorem 3.10), with the same 
poles but different residues: 

W*W=J2x^.^ (3.52a) 

~ -^-1 b* 

W*iX) = -b*^+J2-^, (3.52b) 

^*(A) = ^+EX^> (3.52c) 

2 — 1 

Z*(A)^^/E^- (3-52d) 

Just like in the proof of Theorem 3.10, it follows from Theorem 2.14 that 

^.^^ A, + fij .^^ Aj + Hj 

SO that Z* and Z* are determined by W* and W* . Moreover, there is the 
following connection between the ordinary Weyl functions and their adjoints. 

Theorem 3.15. Assume that K >2. The residues of W and W* satisfy 

aual = ^ -, k=l,...,K. (3.54) 

z— 1 ^ 

Likewise, the residues of W and W* satisfy 

bkbl = ^^ ^' A:=l,...,i^-1. (3.55) 



2n 1 



.=1 ^ ^^ 



(The empty product appearing when K = 2 should be omitted; thus, bibi 
iMi 11^=1(1 + A^i/'^j) *^ ^'**'^ case.) Moreover, 

/K-l s I / K 



'-"-iiT^nn'-'l/inM- p'^^) 
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Proof. We first prove (3.54). From (2.16) we have 53i(-A) = 5'2i(A)532(A) - 
S'22(A)S'3i(A). Evaluation at A = A^: kills S'31, so 

'S'3i(— Afc) = >5'2i(Afc)532(Afc). 

Since the poles of W and W* are simple, the residues are given by Uk — 
-5'2i(Afc)/S'^i(Afc) and al = -S'32(Afc)/S'^i(Afc). Multiplication yields 

* _ 'S'2i(Afe)5'32(Afe) _ 5'3i(— Afc) 

and insertion of the expressions for S'31 and S'31 from (3.39) finishes the job. 
The proof of equation (3.55) is similar. 
As for (3.56), we saw in (3.42) that 

Ooo = 



hx + hx-i 



In the same way, or by using the symmetry transformation (3.62) described in 
the next section, one shows that 

°° lo + h 
highest coefficient of S'31 yields 



Combining S3i(A) == ~2Ani=i(l ^ A/A^) with the expression (3.22) for the 



II ^» ^ 2 II ^^" ) II ^"'^"' 

i=l \m=0 / \a=l / 



and similarly we find by comparing S'3i(A) = ^2AJ^- -^ (1 — ^^/fJ^j) to (3.36) 
that 

K-l ^ /K-l \ /K-l \ 

H Mj = 7^{hK + hK-l) W hm-l {lo + h) [ Y\. 9a+lha ■ 
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(3.57) 

as shown in (4.53), while (3.56) breaks down for the same reason that (3.42) did; 
by (4.49), (4.50) and (4.51), we have instead 

in this case. 
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i=i 


\m=2 




Equation (3.56) follows. 






Remark 3.16. When K = 


= 1, we have 






flifli = 


2 

^ A^ 



Remark 3.17. Theorem 3.15 shows that W and W together determine W* , 
since a|, . . . , a*j^ can be computed from (3.54) if one knows {ak,bk, boo, Ai, fij}. 
But they only almost determine W*; the residues 6^, . . . , b*j^_i can be computed 
from (3.55), but the constant 6^ is not determined! This turns out to be highly 
significant for the inverse spectral problem: the Weyl functions W and W don't 
contain enough information to recover the first weight gi and its position j/i; for 
this we need to know the value of b't^ as well. 

We can now prove the positivity of the residues ai and bj in Theorem 3.10. 
(The notation introduced in this proof will not be used elsewhere, and is omitted 
from the index of notation in Appendix C.) 

Proof of Theorem 3.10, continued. We consider the residues {fli}^]^ first. For 
if = 1 we have S'2i(A) = —giloX and 531(A) = — 2A + gihiloh^^, so that 



W{X) 



521(A) ^ 
'531(A) A 



1 

hih 



gihiloh 



hence ai — 



hih 



> 0. We now proceed by induction on K. Suppose that the 



residues Ui are positive when K = m — 1, and consider the case K — m > 2. 
Because of (3.54), no a^ can ever be zero as long as all masses are positive, and 
therefore it is sufficient to verify that all ai are positive when the last pair of 
masses are given by gm — h,n — £ with e > small; since the residues depend 
continuously on the masses, they will keep their signs as gm and hm are allowed 
to vary arbitrarily over all positive values. From (3.9) we get 



'5ii(A,£ 

521 (A, £) I =i2m(A) 
,53i(A,£ 



L 



2m- 1 



(A) 



J2m-2 



(A). 



Li(A) 



io(A) 



where we consider all positions and all masses except gm = hm = e as fixed, and 
treat the 5y(A,e) as polynomials in two variables. The spectral data defined 
by these polynomials will then of course also be considered as functions of e: 
{Xi{e),ai{e)}^i. (As we will soon see, the largest eigenvalue A„j(e) has a pole 
of order 2 at e = 0, while the other eigenvalues are analytic functions of e.) The 
first four factors in the product above are 



i2m(A) 



Lo 



l(A) 



1 


-{hm + ^2m-l)A 



£ 

1 

— shmX 



e 
e 
1 ~ £ hmXj 



We denote the product of the remaining factors by (sii(A), S2i(A), S3i(A))-^; these 
polynomials have the same form as 5ii, 52i and 523 (see Corollary 3.4), but 
with m — 1 instead of m, so their degrees are one step lower, and they only 



depend on {gk,hk}™^i and {h} 
thus have 



2m- 
k=0 



not on hm-i, hm and g„ 



We 



'5ii(A,£) 

52l(A,£) 
i53i(A,£) . 



<h 



hm-l)X 

%i(A,0)\ /O 

52i(A,0) + 

.53i(A,0)/ VO 




(3.59) 
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The polynomials ^^^(A, 0) define the spectral data for the case K — m — 1 
(since the final pair of masses is absent when e = 0); in particular we know 
from Theorem 3.8 that 5'3i(A, 0) has a zero at A = 0, and that the other m — 1 
zeros are positive and simple. If A = A^ 7^ is one of these other zeros, then 
at the point (A,e) = (Ai,0) we therefore have S31 = and 9S'3i/i9A 7^ 0, so 
by the Implicit Function Theorem there is an analytic function A^ (e) , defined 
around e = 0, such that Ai(0) = A^ and S'3i(Ai(e),e) = 0. It follows that for 
i = 1, . . . ,m — I, the residue 

ai{e) = res W{\,e) — 



depends analytically on e too, and it is therefore positive for small e > 0, since it 
is positive for e = by the induction hypothesis. This settles part of our claim. 
It remains to show that the last residue am(e) is positive. As a first step, we 
show that \m{£) has a pole of order 2 at e = 0. For convenience, let 

/(A,e) = y^ — ; 

this is a polynomial of degree m in A, and Am(e) is the largest root of the 
equation f{X,e) = 0. From (3.59) we have 

/(A,£) = /(A,0)-?2m(es2i(A) + e's3i(A)). 

Using Corollary 3.4, we see that the leading terms of /(A,0) = S'3i(A,0)/A and 
^2mS3i(A) are (-l)™CiA™-i and (-1)™C2A™, respectively, with 

(m-2 \ /m-1 \ 

Y[ ^2r (^2m-2 + hm-l + hm) J^ QaK > 0, 

Cm \ /r?? — 1 \ 

n '2r J I n 9aha I > 0. 

(The precise form of these constants is not very important, only their positivity.) 
Moreover, S2i(A) has degree m — \. Thus 

/(A, e) = /(A, 0) - hm (£S2i(A) + e^s^^iX)) 

= (-l)™+iC2e2A™+p(A,e), 

with a polynomial p(A, e) of degree m— 1 in A. Since p{\, 0) = /(A, 0) has leading 
term (-l)"CiA"-i, we see that 

e2"-2p(^g-2^^) ^ (-1)™Cik'"^i + (terms containing e). 

Hence, the equation /(A,e) = 0, of which A„j(£) is the largest root, can be 
written in terms of the new variable k — \e^ as 

0=(-l)™+l£2™-2j(^^£) 

= Cae^^A" + e2'"-2(_i)™+ip(A, e) 
= C2K" + e2™-2(_i)"+ip(«;e-2^ e) 
= C2«;"-Cik"-i + £?(«:,£), 
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for some two- variable polynomial q{K,,e). As before, the Implicit Function 

Theorem shows that this equation has an analytic solution ^(e) with k(0) = 

C1/C2, which corresponds to a meromorphic zero of /(A, e) with a pole of order 2, 

as claimed: 

, , ,^^{e) _ Ci/C2 + 0{e) 



Finally, the corresponding residue is 



OmCs) = res W(\,e) 

A=A„(e) 



5'2i(A™(e),e) 






(A„(e),e) 



The derivative of the polynomial 531 at its largest zero has the same sign as 
the leading term of 531, namely (—1)™+^. As for the sign of ^21, we have from 
(3.59) that 

^2i(A,£)-52i(A,0)+es3i(A), 

where 5'2i(A, 0) and S3i(A) have degrees m— 1 and m, respectively. When this 
is evaluated at A = Am(e) ~ §^e~^, the two terms on the right-hand side are 



of order e^™ ^ and e^™ ^ , respectively, so the dominant behavior as e 
comes from the leading term of S3i(A): 

C2 



0+ 



^2l(A™(£),£)~ £(-!)' 



h 



C1/C2 



In particular, the sign of S'2i(A„i(£),£) is (—1)™, and it follows that am{£) > 0, 
which is what we wanted to show. This concludes the proof of positivity for the 
residues at. 

The proof for the residues {bj}jSi^ is similar. In the base case K — 1 there 
is nothing to show. Assume that they are positive for K = m — 1, and consider 
the case K — m> 2. We have from (3.23) 



^ll(A,£)' 

^2l(A,£) 

.53l(A,£), 



i2m(A) 



L. 



l(A) 



.(A)' 



Li(A) 




Splitting off the first four factors 



L2m{X) 



we obtain 



L 



2m- 1 



(A) 



1 

-{hm + hm-l)^ 



I — S hm-l^ 
"e(^2m + ^2m-l)A 1, 



'^ii(A,£)\ /5ii(A,0)\ / e 0\ /sn(A)^ 

^2l(A,£) = ^2l(A,0) + -£/2m-lA -£'/2m-lA £ S2i(A) 

^S'3l(A,£)/ \S'3l(A,0)/ \ -£(;2m+^2m-l)A O/ Vs3l(A), 

(3.60) 
where the degrees on the left-hand side are (m — l,m, ?n), while both 3x1 
matrices appearing on the right-hand side have degrees (m — 2, to — 1, m — 1) 
(cf. Corollary 3.7). The eigenvalues {fJ'j{£)}Tj~i^ are the zeros of the polynomial 



/(A,£) 



53l(A,£) 53l(A,0) 



A 



A 



^{hm + '2t?j--i)s2i(A). 
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As above, it follows easily that {Mi(£)}7Li^ ^^'^ analytic, and that the corre- 
sponding residues {bj{e)}"lr^ are positive. The largest zero ^m_i(£) has a pole 
of order 1 at e = 0, as we now show. By Corollary 3.7, the leading terms of 
^3i(A,0) and {hm + l2rn~i)s2i{X) are {-l)"^'^CiX"^'^ and {-l)"'-^C2X'^~\ 
respectively, with some positive constants Ci and C2. (For the record, these 
constants are 

/ 2m \ /m-2 \ /m-2 \ 



\ r^2 / \ a^l 

771—1 \ /m—2 



C2 = {hm + hTTi-l) Y\. ^^r-l (^0 + ^1) h,n-l Y\. da+lK > 0. 

Special case: Ci = 2 if to = 2. The empty product rir=2 '2r-i is omitted in Ci 
when m — 3 and in C2 when to, = 2.) Hence, 

7(A,e) = /(A,0) - eihrn + l2ra-l)s2l{X) 

= (-l)"C2eA"-i+p(A,e), 

with a polynomial p{\, e) of degree m — 2 in A, such that p(A, 0) = /(A, 0) has 
leading term (-l)™-iC'iA™-^ so that 

£'"-2p(Ke-i,e) = (-l)"-!^!^""^ + (terms containing e). 

The equation /(A,£) = 0, of which /^m„i(e) is the largest root, can therefore be 
written in terms of the new variable k — Xe as 

0=(-l)™£™-27(A,£) 

= C2e""^A'"-i + e"-2(_;L)™p(A, e) 

= C2K™ - CiK™-! + e q{K, e), 

for some two- variable polynomial q{K,e). The Implicit Function Theorem gives 
an analytic function K(e) with k(0) — C1/C2, and, as claimed, 

,, K{e) Ci/C2 + 0{e) 
Hj-i{e) = ^ = . 

The corresponding residue is 

5'2i(Mm-i(e),e) 



6m_i(£) = res W{X,e) = — 

The leading term of S'31 determines the sign of the derivative dS^i/dX at the 
largest zero, namely (—1)™. From (3.60), 

S'2i(A,e) = 52i(A,0) -d2m-iAsii(A) - e^/2m_iAs2i(A) + es3i(A), 

and when evaluating this at A = /i„j_i(e) ~ W-e^^, the last three terms on the 

C2 

right-hand side are of order e^~™, so the contribution of order e^~™ from the 
first term S'2i(A, 0) is the dominant one as £ -^ 0+, and it has the sign (—1)™^-'^. 
It follows that 6m_i(e) > 0, and the proof is complete. D 
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3.6 Symmetry 

For solutions of the differential equation (2.2a), ^ = yl(A)$, the transition 
matrix S{X) propagates initial values at the left endpoint to final values at the 
right endpoint: 

*(+f) = 5(A)$(-f). 

The transition matrix depends of course not only on A but also on g(y) and 
h(y), which in our discrete setup means the point masses gj and hj interlacingly 
positioned at the sites yk with Ik — yk+i — Uk', let us write 



S{\) = L2k{\) 



hK 




L 



2K-1 



(A) 



■ 0' 

.9K. 



— S{X;lo, . . . ,l2K]9i,hi 



L2K-2W ■ 
9k, hK) 



ii(A) 



io(A) 



to indicate this. 

For the adjoint equation (2.22a), ^ = y4(— A)fJ, we saw in Proposition 2.f2 

that the matrix S'(— A)^^ = J S{X)^ J propagates values in the opposite direction, 
from initial values at the right endpoint to final values at the left endpoint. If 
we denote this matrix by S*{\), we thus have 

n{~\) = S*{\)Vi{+l). 

When going from right to left, one encounters the point masses in the opposite 
order compared to when going from left to right, and the following theorem 
shows that the solution ^{y) reacts just like $(?/) does when encountering a 
mass, except for a difference in sign. 



Theorem 3.18. The adjoint transition matrix is given by 

S*{X) — S{X; I2K, ■ ■ ■ Jo', —hK, —gii, ■ ■ ■ ,—hi, — ffi)- 
(And similarly with tildes for the twin problems.) 



(3.61) 



Proof. Use J ^ J^ ^ J-^ together with JLk{\fJ = Lk(\) and J[^J\^ J = [l^] 



to obtain 

S*{\)^JS{\fj 

^J\L2k{\) 



Hk 




JLoiXfj 
io(A) 



-.91 
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-Hk 



io(A) J 



J)...[J 

L2k{X) 



Hk 




j]ijL2K{Xfj 



D 



Remark 3.19. The adjoint Weyl functions W* and Z* are defined from the 
first column in S* = JS^J, 

{Sii,Sli,Sli) = (533, —5*32, 5*31) , 
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in almost the same way as W and Z are defined from the first column in S, but 
there is a slight sign difference in IF* since we have defined all Weyl functions 
so that they will have positive residues: 

W = -821/831, Z = -811/831, 

but 

W* = —832/831 = +S21/831, Z* = —833/S31 = —8I1/831. 

As a consequence, we see for example that if 

flfe — Filo, ■ ■ ■ , hK',gi, hi, . . . ,gK, hx) 

indicates how the residue a^ in W depends on the configuration of the masses, 
then 

^'^fe — F{hK, ■ ■ ■ ,la', —h-K, —gx, ■ ■ ■ , ^hi, —gi), 
with the same function F , will determine the corresponding residue in W* . 

Remark 3.20. In Section 4.2 we will use Theorem 3.12 and Proposition 3.14 
to derive formulas for recovering the weights hj and their positions j/2j from the 
Weyl functions W and W . Because of the symmetry properties described here, 
the same formulas can then be used to recover the weights gj and their positions 
2/2J-1 from the adjoint Weyl functions W* and W* , by substituting 



(3.62) 





Ci i-^ -a*. 


h ^ hK-k, 




b,^-b*. 


gj ^ -hx+i-j; 




boo ^ -bla. 


hj ^ -gK+i-j- 


Note that 1 — j/„j 


-Efcf™^Ustobe 


replaced by 22k=, 


1 + y2K+l-m- 







J \-^2K—m 

hK-k — 2^s=o ' 



4 The inverse spectral problem 

To summarize what we have seen so far, the Weyl functions Vt^(A), Z{\), W{X), 
Z(X) encode much of the information about our twin spectral problems. In 
particular, in the discrete interlacing case with positive weights, the Weyl 
functions are rational functions in the spectral variable A, with poles at the 
(positive and simple) eigenvalues of the spectral problems, and the functions 
Z and Z are completely determined by W and W (which in turn are of course 
determined by the given discrete measures m and n that define the whole setup). 
The measures depend on the AK parameters 

Xi,X2,- ■ ■ ,X2K-1,X2K, mi, 1713, . . . , m2K-l, n2,nA, ■ ■ ■ ,n2K 

(or equivalently {j/fc, 520-1, /i2a}), while the Weyl function W depends on the 
2K parameters 

Ai, . . . ,\k-, ai,. . . ,aK, 

and its twin W on the 2K — 1 parameters 

fJ-i, ■ ■ ■ ,fJ-K-i, bi,. . . ,bK-i, boo- 
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To get an inverse spectral problem where the number of spectral data matches 
the number of parameters to reconstruct, we therefore need to supplement W 
and W by one extra piece of information, and a suitable choice turns out to be 
the coefficient b^ defined by (3.52b). We will show in this section how to recover 
the discrete interlacing measures m and n (or, equivalently, their counterparts g 
and h on the finite interval) from this set of spectral data {Aj, aj, /Zj, bj, b^o, b*^} 
that they give rise to. Moreover, we will show that the necessary constraints 
(0 < Ai < • • • < Xk, < /ii < • • • < ^JLK-i, and all a^, bj, boo, b'^ positive) are 
also sufficient for such a set of numbers to be the spectral data of a unique pair 
of interlacing discrete measures m and n. 

4.1 Approximation problem 

As we mentioned in Section 3.4, the properties in Theorem 3.12 are enough to 
determine the polynomials Q, P, R uniquely, and this fact will be proved here. 

Theorem 4.1. Let boo be a positive constant. Let a and (3 be compactly supported 
measures on the positive real axis, with moments 



ak= x^da{x), 13k = / y^dPiy), (4.1) 

and bimoments (with respect to the Cauchy kernel jn— j 

hm^ // -^da{x)d(3{y). (4.2) 

JJ x + y 

Define W , W , Z , Z by the formulas (3.44) (repeated here for convenience) : 

(3.44a) 

(3.44b) 

booW{X), (3.44c) 

(3.44d) 

Fix a positive integer j . (If a and f3 are supported at infinitely many points, then 
j can be arbitary. In the discrete case with a = X]i=i '^i^Xi o-nd (3 — X]i=i biS^., 
we restrict j to the interval 1 < j < K .) 

Then there are unique polynomials Q{X) — Qj{X), P{X) = Pj{X), R{X) = 
Rj{X) satisfying the conditions of Theorem 3.12 (also repeated here for conve- 
nience): 

degQ = j, degP = j-l, degR = j-l, (3.46) 

g(0) = 0, P(0) = 1, (3.47) 



/C [A) = 


J X-x' 


^(A) = 


f dfi{y) , 
J A-y 


Z(A) = 


1 ff da{x)dp{y) 
2X^JJ {X-x)ix + y) 


Z(A) = 


1 ff da{x)dl3{y) 
2A J J {x + y){X-y) 
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and, as A — > oo, 



W{X)Q{X) - P{\) = O 
Z{X)Q{\) - R{X) = O 



i?(A) + P{X)W{-X) + Q(A)Z(-A) = O 



XJ 



These polynomials are given by 

Q(A) = Ap(A), 

^(,) /Q(A)-Q(.) 

J X-x 

Q(A) - Q{x) 



R{X) = 



{X-x){x + y) 



da{x)dP{y) + \p{X) + booP{X), 



where 



det 



p(A) = 



det 



/ 1 ho 

X /20 

( OtQ Iw 

Oil ha 

\"i-i ho 



hj-2 



hj-2/ 

h,-2\ 

h.j-2 

h,J^2j 



(If j = I, equation (4.3d) should be read as p(A) = l/aQ.) 
In particular, we have (using notation from Section A. 3) 



det 



Q'(0)=p(0) = 



'20 



ijO 



L 



2J-2 



h.j-2, 



det 



/ ao ho 
ai ho 

\"j-i ho 



h.j-2\ 

h,j-2 



^JJ-2/ 



-t2U 
-tOI 



(3.48a) 
(3.48b) 
(3.48c) 

(4.3a) 
(4.3b) 

(4.3c) 



(4.3d) 



(4.4) 
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(to be read as Q'(0) = l/ao if j = I), and 
p{x) 



i?(0) 



da{x)d(3{y) + ip(0) + b, 
x + y ^ 

I loo +2 -'^lO ■ 



det 



'10 



'20 



V ^J-1,0 IjQ 



^2,] -2 



IjJ-^J 



det 



/ ao ho 
ai ho 



h,,-2\ 
h,3~2 



6oo = 



/C," 



(4.5) 



'tOI 



+ 6o 



\"j-i ^jo ■ • • hi-"^) 
(to be read as i?(0) = {ho + \)/<^o + ^oo if j = Ij- 

Proof. A bit of notation first: define projection operators acting on (formal or 
convergent) Laurent series /(A) — J^kez CfeA'" as follows: 

Il>of = Y,CkX'', ^>o/-^CfeA^ n<o/ = ^CfcA'=. (4.6) 

k>0 k>0 fc<0 

Note that we can expand Vl^(A) in a Laurent series with negative powers, 



W{\) 



da(x) 
X — X 



E&M^^-T^. 



k>0 



fc>0 



and similarly for the other Weyl functions. 

We see at once that the conditions (3.48a) and (3.48b) determine the polyno- 
mials P and R uniquely, by projection on nonnegative powers, if the polynomial 
Q is known: 

P = Il>o[QW], i? = n>o[gZ]. (4.7) 

Inserting this into (3.48c) gives 



A^ 



n>o[QZ](A) + n>o[QW^](A) W{-X) + Q(A)Z(-A) = O 

Writing n>o = id — n<o produces 

Q(A)Z(A) + Q{X)WiX)Wi-X) + Q(A)Z(-A) 

- n<o[QZ](A) - n<o[QW^](A) W{-X) = O 



Xi 



where the first three terms cancel thanks to the identity Z{X) + W{X)W{—X) + 
Z{—X) = which follows from the definitions (3.44) by a short calculation (cf. 
also (2.19)). This leaves 



n<o[QZ](A) + n<o[OH^](A) W{-X) = O 



X3 



(4., 
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Next, note that 

Q{X)W{X) ^ Q{X) /^= f ^('j - ^("^ da{.) + f^dai.), 
J \~x J X — x J X — x 

where the first term is a polynomial in A (since Q{X) — Q{x) vanishes when 
A = a; and therefore contains A — x as a factor), and the second term is O (1/A) 
as A — > oo. Thus, the first and second term are n>o[QW^] and n<o[QW^]j 
respectively, which gives on the one hand the claimed integral representation 
for P, 

P{X)=Il^,[QW]{X)^ j^^-^^da{x), (4.9) 

and on the other hand, multiplying the negative projection by W{—X), 

n<.[<3»1(A) m-x) = (/ ^ M^)) (/ ^ - "o 

Similarly, 



^^^^ -da{x)d(3{y) + b^Q{X)W{X), 



{X-x){x + y) 

where the first term is a polynomial in A since we require (5(0) = 0, likewise the 
second term is a polynomial (by the same argument as above), and the third 
term is C'(l/A) as A— >oo. Thus, we obtain from the first two terms, together 
with the contribution to nonnegative powers from the fourth term, the claimed 
integral representation for R, 

R{X) ^ n>o[QZ](A) = 9i^ + JJ ^f^l'^^^^^l^ da{x)d/3{y) + b^P{X), 

(4.11) 
and from the third term, together with the contribution to negative powers from 
the fourth term, 

n<o[QZ](A) = // ,. ^\f^ , da{x)d/3{y) + b^ f ^ da{x). (4.12) 
J J {X-x)(x + y) J X-x 

Inserting (4.10) and (4.12) into (4.8) gives 



which simplifies to 



^^""^ da{x)dp{y)^o(^]. (4.13) 



{x + y){X + y) '^^'" \Xi 
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Since Q(0) = 0, we write Q{x) = xp{x), where p is a polynomial of degree j — 1. 
Upon expanding 1/(A + y) = X]a;>o y'"^ '■'^^^■'' ^^'^ condition (4.13) takes the 
form 

pM^^ da{x)dp{y)^0, 0<k<j~2. (4.14) 

x + y 

This imposes j — I linear equations for the j coefficients in p{x) = pq + pix + 
\-pj-ix^^^: 



(po,---,Pj-i) 



fho 

I 20 



V: 



jo 



^2 J -2 



^J,i-2/ 



(o,...,o). 



Adding an extra column, 



iPo,---,Pj-i) 



fho 

1 20 



V: 



jo 



h,j-2 h.j-i\ 

l2,j-2 h.j-l 



-^ij-2 Ij.3-1/ 



(0,...,0,*), 



we see that the row vector (po, • • • tPj-i) is proportional to the last row of the 
inverse of the bimoment matrix in question (which is invertible by the assumption 
about infinitely many points of support, or by the restriction on j in the discrete 
case). Hence, by Cramer's rule. 



p{x) — Cdet 



(ho 
ho 



h.3-2 1 \ 

h.'j-2 X 



I 



0,3 



-2 X^-'J 



The constant C is determined by the remaining normalization condition P(0) = 1; 
from (4.9) we get 



1 = P(0) = j 



Q(0) - Q{x) 



Cdet 



(ho 
ho 



Q-x 



da{x) — / p{x) da{x) 



j-2 

^2 J -2 



Jlda{x) \ 
J xda{x) 



Cdet 



(ho 
ho 



\ho 



Ij,]-2 J x^ ^da{x)J 



h,j-2 ao \ 
h,j-2 oil 

4j-2 aj-i) 



Finally, the expression (4.4) for Q'{Q) follows at once upon setting A = in 
Q'{\) = p(A) + Ap'(A) and using the determinantal expression (4.3d) for p(A); 
the last term in (4.4) represents an evaluation of the determinants in terms of 
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certain integrals (which will be sums when the measures a and /3 are discrete). 
This is explained in the appendix; see Section A. 3, in particular equations (A. 18) 
and (A. 21). The expression (4.5) for i?(0) is also immediate from the formula 
(4.3c) for R{X), since (g(0) - Q{x))/{0 - x) = Q{x)/x = p{x). (The symbol ICj 
is just notation for the determinant in the numerator; it doesn't seem to have 
a simple direct integral representation, but we will mainly be interested in the 
difference between Rj{0) and i?j+i(0), which equation (A. 27) takes care of.) D 

Remark 4.2. The polynomial p{x) in Theorem 4.1 is proportional to pj^i{x), 
where {pn{x) , qn{y)} n>o are the normalized Cauchy biorthogonal polynomials 
with respect to the measures xda{x) and dj3(y). This can be seen either directly 
from the biorthogonality condition (4.14), or by comparing the numerators in 
the formula (4.3d) for p and the formula (A. 4) (with Ia+i,b instead of lab) for Pn- 

4.2 Recovery formulas for the weights and their positions 

Most of the work is now done, and we can at last state the solution to the inverse 
problem of recovering the weights gj and hj and their positions yk from the 
spectral data encoded in the Weyl functions. The answer will be given in terms 
of the integrals J^m defined by equation (A. 15) in Section A. 3 in the appendix. 
Since a and /3 are discrete measures here, these integrals are in fact sums; see 
(A.31) in Section A.4. 

Theorem 4.3. The weights and positions of the even-numbered point masses 
are given by the formulas 

hj,^^J^^±l + b^, (4.15) 

ao 

{l-y2K)hK^—, (4.16) 



ao 



and, /or j = 2, . . . , K, 





nrW ( 'TOO 1 1 -Til \ 


-tOI -tOI 

-Til -tIQ 
-Jj-ld-2Jj-lJ-l 

-tOI -tOI 
^j-l,j-2^j,j-l 



(4.17) 
(4.18) 

Proof. We use Proposition 3.14 together with (4.4) and (4.5) from Theorem 4.1. 
The rightmost mass is special; we have 

hK = i?i(0) - RoiO) = (l^^^±A + bJ)-0 



1 
ao 



ao 
and 

{I - y2K)hK = Q[{0) ~ Q'oiO) 

For the other masses we get 

hK+i-j = Rj{0) - Rj-m = 



^] , u \ i ^1- 



-tOI I '^oo I I „oi 
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which equals (4.17) according to (A. 27), and 

-t20 -t20 



(1 - y2iK+i-j))hK+i-j - g -(0) - q;_i(o) 



-tOI -tUI ' 

which equals (4.18) according to (A. 24). D 

The symmetry described in Remark 3.20 immediately provides formulas 
for the odd-numbered point masses. We let {J'*)nm denote the integral J^^ 
evaluated using the measures 

K K-l 

a*=Y,<h. and /?* = ^ 6*5^^. (4.19) 

in place of a and /3, and similarly for the moments a*. = {J*)io and /3* = {J*)'q\, 
and the Cauchy bimoments I*^ = (i/*)^f . Then the symmetry transformation 
(3.62) also entails the substitution 

XL ^ {-^T^'^iJlnr.. (4.20) 

(including as special cases a^ i— > — a^, Pk H' — /3^, and /afc i— > /*j^). 

Corollary 4.4. T/ie weights and positions of the odd-numbered point masses 
are given by the formulas 

91 -'^ + b*^, (4.21) 



a 







{l + yi)gi = ^, (4.22) 

an 



^0 

and, for j ^2,...,K, 

/'-T*\io /'('-7*^oo I i/-T*\ii ^ 

w Ji-lj-2W lj,j-l 

( -7* "ill (" -7*';io 

(1 + y..-i)5. - ^f^b■^-^^fi^^^-^ (4.24) 

Corollary 4.5. The corresponding weights and positions on the real line are 
given by the formulas 

X2K = ^ In 2(/oo + booao), (4.25) 



1 / Joo + booOlo ,. „„. 

n2K = — V ^ , (4.26 

ao V 2 



.xi = -iln2(/*o + &:o"S), (4.27) 



,„, __ LJIk±^^ „..s, 
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and, for j = 2,..., K, 



X2(K+1- 



J) 2 



/ 2 7-00 






n2{K+l-j) - -7m 



qW / 'TOO -711 



-tOI -tOI 

>-'.7-l,l"-2>-' 7,7-1 



(4.29) 
(4.30) 



a;2j-i = -5 In 



'2(J*^oo 

(:^"i 



jj-i 



j-i,i-2. 



TO2J-1 



('^*)j-lj-l , / (^*)ij-l(^*)j-lj-2 



w ;7-i,7"-2W Jjj-i 



(4.31) 



(4.32) 



/n terms of non-starred quantities (together with b*^), the odd-numbered variables 
take the form 



( 2J^P" ' 
X2{K+i-])-i = s^n -xii ■'■' — 
V'-'i- 1,7-1 



m2(K+l-3)-l - ^10 ^10 



-rOl / -Til -TOO 

'-'.7,7-1 . >-'7-l,7-W.7.^ 



'J-lJ-l'-JJ 



(4.33) 
(4.34) 



for j — 1, . . . ,K — 1, and 
( 



2^"o 



K.K-l 



-Til 



K-2 



2^-Tio , 



M/L 



-Tio \ 2 

•-'if-i.^r-i \ ^ 



-TOO 



K.K-l -Til 

>-'K-l,i<'-2 



^^^10 



(4.35) 



(4.36) 



where L = Y\i=i ^i '^'^d M = Jli^i Mj- 



22;fc 



les 



Proof. Since j/fc = tanhx^ = (e *= — l)/(e *= + 1), we have 

exp(2xfe) = = 1. 

1 - y/c 1 - yfc 

Moreover, hj — 2n2j coshx2j — n2j{e'^^^^ + l)e^^^^ = n2j ( -^_^ J e^^^J inipl: 
that 

2n2j exp(-X2j) = (1 - 2;2j)^i- 

Now it is just a matter of plugging in the formulas from Theorem 4.3 and solving 
for even- numbered X2j and n2j. For example: 

2exp(2a;2(K+l_j)) = = ry 

l-J/2(K+i-i) 2 {l~y2{K+i-j)jhK+i-j 2 

-Tio ( -TOO I 1 -Til ^ -1 -TOO 

_ ^j-ij-iWi,j-i + 2^j-ij-2J -1- _ ^j,i-i 



-Til -TlO 

>^J-l,i-2>^i-lj-l 



2 ^7^1 



j-i,i-2 
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The odd- numbered X2j-i and m2j-i are dealt with similarly, using the formulas 
from Corollary 4.4 together with 

exp -2xfe = — = — 1 

1 + 2/fc 1 + J/fc 

and 

2to2j_i exp{x2j-i) = (1 + y2j-i)gj- 

In order to translate starred to non-starred, we use Lemma A. 3 (with A = K 
and B = K -1): 

j^2j-U-l)+0-ljyp(j-l)-j+Q-l -T.1-04-0 

o2?"-l -rOO ' 

^ ^K,K-\ 

and similarly for the other {^J*Ynm occurring in the formulas for X2j-\ and 
TO2j_i. All the factors L, M and J'^}^_x cancel in the quotients, except in the 
formulas for x\ and m\ where we have 



'K,K-\ 

and 



«s = (:^*)i§ = {r)\ 



10 — W 710 — 2/7°° 



n 



Remark 4.6. A more compact way of writing the solution is to state the 
formulas in terms of the following quantities (where r — K + I ~ j throughout) : 

i exp 2x2K = ^0 + booao, 
^exp2a;2r = ^i"~^ ■- i = 2,...,X, 

jOO 

i exp 23^2.-1= „"" , j = l,...,X-l, (4-37) 



-TOO 
1 o '^' 

i expzxi 



K,K-1 



2b^L 
M 



-Til I OO -TlO 

^K-l.K-2 ' 7i,,r '-'Ar-l.K-l 



and 



2n2_R-exp(-a;2i^) = — , 
ao 

-Til -TlO 

2n2rey.Y>{-X2r) = ^hfoT -^ — ' J = 2,..., if, 

7.11 TPl 

2m2r-iexp(-a;2r-i) = 'Tio'^io '^"^~^ ' J = l,---,^-l, 

rll 



(4.38) 



^^ ' ^K-l,K-2 

TlO 

'K-l,if-l 



2mi exp(-a;i) = .- ^^o ' + 26^ 
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We now know that the set of spectral data computed from the interlacing 
discrete measures m and n allows us to reconstruct these measures uniquely, and 
we also know (Theorem 3.8, Theorem 3.10, equation (3.56)) that the eigenvalues 
are positive and simple and that the residues are positive (provided that the point 
masses in m and n are positive). Next we will show that there are no further 
constraints on the spectral data, i.e., any set of such numbers are the spectral 
data of a unique pair of interlacing discrete measures. It will be convenient to 
introduce a bit of terminology first. 

Definition 4.7. Let V C R"*^ (the "pure peakon sector") be the set of tuples 
p = (xi,... ,X2K;mi,n2,... ,m2K-i,n2K) 

satisfying 

xi < ■■■ < X2K: all m2a-i > 0, all n2a > 0, 

and let TZ C R''^ (the "set of admissible spectral data") be the set of tuples 

r = (Ai, . . . , A/f;/ii, . . . ,^K-i]ai, . . . ,a_R-; 6i, . . . ,hK-i;boo,h*^) 

satisfying 

0<Ai<---<Ak, < /xi < • • • < /ik_i, all aj,6j,6oo,^|^ > 0. 

The forward spectral map taking a point p S 7^, representing a pair of 
interlacing discrete measures 

K K 

m^2'^m2a-i5x2a-i and n = 2 y^ n2a 6^2^ , 

a—l a—1 

to the corresponding spectral data r G TZ (as described in Section 3) will be 
denoted by 

S:V^n. (4.39) 

The formulas in Corollary 4.5 (or Remark 4.6) define a function 

T-.n^R'^^ (4.40) 

which we will call the inverse spectral map. 

Theorem 4.8. For K > 2, the function S maps V bijectively onto TZ, and 
T '■ TZ ^ T' is the inverse map. (See Section 4-3 for the case K — \.) 

Proof. To begin with, T maps TZ into T'] this is the content of Lemma 4.9 below. 
By Corollary 4.5, T o S — id-p. Thus 5 is a homeomorphism onto its range. It 
remains to show that the range of S is all of TZ and that S oT = id^. For this, 
it is most convenient to use the alternative description of the forward spectral 
map given in Appendix B, where the spectral data are defined directly in terms 
of {xk, m2o-i, 'n'2a} without going via the transformation to the finite interval 
[—1,1]. (At the beginning of Appendix B there is a summary comparing the 
two descriptions.) The eigenvalues Ai and /ij, as well as the residues a^, bj 
and boo, are all uniquely determined by certain polynomials ^(A), A(X), B(X) 
and -B(A) with the property that their coefficients are polynomials in the variables 
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{m2a-ie=t^^"-Sn2ae=^^^°}; see (B.ll), (B.12), (B.19) and (B.20). Let us write 
A{X;p) (etc.) to indicate this dependence of the coefficients on the masses and 
positions. It is clear from the symmetry of the problem that 6^ could also 
be definied similarly (although we have chosen not to name and write out the 
corresponding polynomials, instead using (B.6) as the definition). 

Now, since Remark 4.6 exhibits the variables {m2a-i 6^^^"-^ n2a e*^^°} as 
rational functions of the spectral variables r, the coefhcients in the polynomial 
A{X; T(r)) are also rational functions of r. Since T o S = id-p, we know that 
these coefficients agree with the coefficients of ni=i(l ~ A/Afc) (see (B.3)) for 
each r in the range of S (which is an open set in TZ since 5 is a honieomorphism). 
Hence A{X; T(r)) = nj=i(l ~ A/Afc) identically as a rational function of r, and 
in particular this identity holds for any r G 7?.. The same argument works for 
the other polynomials involved in defining the spectral variables, and therefore 
5o T = id^, as desired. D 

Lemma 4.9. The function T maps TZ into V , i.e., the formulas in Corollary 4.5 
give positive masses m2a-i > and n2a > 0, and ordered positions xi < ■ ■ ■ < 
X2K, for any spectral data in the admissible set TZ. 

Proof. Positivity of m2j^i and n2j is obvious. To show that the positions Xk 
are ordered, we will use the formulas (4.37) for Qk = 2 exp2a;fc and show that 
qi < ■■ ■ < q2K. 

The outermost intervals present no problems, since 



q2K - q2K-l = V^OO -t- "oo"0j - 

and 



q2K - q2K-l = {IqO + bocao) - — rjY = ^00^0 > 

'-'00 



26* L 

1 1 ^ ^K-UK-2 + -^ Jk-1,K-1 JkU,K-2 ^ ^b*^L Jk^K-l ^ „ 

„ „ -TOO -TOO ]\J -TOO ^ ^' 



As for the other distances, the differences 



and 



q2r - 92r-l — ^^rjj ^TTT ' J = 2, . . . , i^ - 1, 

-TOO -T( 

-Til -Tl 

•^j-ij-i "^ij-i 



fer-l - q2r-2 = ^n " ^n '^ , j = 1, • • • , -f^ " 1, 



are all positive according to (the fairly technical) Lemma A. 4 in Section A. 4. D 

Example 4.10. Let us explicitly write out the solution formulas for the inverse 
problem in the case K — 2,hy expanding the sums J,j4 (including /qo = Jii) 
as explained in Section A. 4; recall that the lower indices n and m give the 
number of factors a^ and bj in each term, and also determine the dimensions of 
the accompanying Vandermonde-like factors ^/j (see (A. 33)), while the upper 
indices r and s are the powers to which the additional factors Ai and /ij appear. 
The spectral data are 

Ai , A2 , /ii , ai , a2 , 6i , 6oo , &00 1 
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and we want to recover 

In terms of the quantities from Remark 4.6, we get 

ke " = /oo + Oooao = T — , + \ — , 1- Ooo(ai + 02), 

TOO 






(A1-A2)' 

aia20i 



:^"_ (Ai + Mi)(A2 + m) 
10 



i/io AiOi + A2a2 

-TOO 

1^2x1 _ >^21 



;e 



M 

(Ai - A2) 



-Til I 00 -TlO 

'-'10 + n,f >^11 



2 



(Ai + /xi)(A2 + ^i) 



26;^AiA2 / Aiai5i , A2a26i 
Aiai + A2a2 + 



Ml VAi+/ii A2 + A*] 



and 



1 



ao ai + a2 
2m e^"^^ — '-^oo '"^^o — ^ ' 'Jio _ "i + °2 



^Ai + Ml >^2+fJ-l, 

Example 4.11. Similarly, in the case K = 3 the spectral data are 

and we want to recover 

xi,X2,X3, 2:4, xs, X6, mi, n2, ma, ^4, ms, ng. 



(4.41) 



JiiJio Jli ■ 1 ^i«i^i ^ A2a26i ' 
Ai + /ii A2 + A^i 
-, , , . ( Aiflifei A2a25i A 

-tOI 701 / n2 ' 

^10-^21 /il(Ai-A2) 

("^+"^^Ai+^0(A2 + Mi)"^"^'^ 

2mie-- - ^^^^ + 26* - Mi(Ai«i + A2a2) 

A1A2 
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The solution is 



3 2 , 






3 2 

e 






2 

tynn ., - 



<^oo i^i ,^1 '^i + Mj 



00 , ,rr„ „„ ^ (-^ii + Mj)(-^J2 + A*j) 



^21 _ -f =12,13,23 i 
'10 



lg2a;4 _ _ 

2 aiiai20ib2 

1 2.3 _ ^ _ 7=12,13,23 n.g/n,-^i(A.+A^j) 

2*^ ~ -rll - 3 2, 



4 = 1 J = 

'2(\. _ \^\2{\^ _ \^\2( , ,,^^2 

0102036162 

^3^2° ntin=i(A.+A^,) 



(Al - A2)2(Ai - A3)2(A2 - A3)2(^i - M2)' 



1 „2a:2 _ '^32 ^ 

2^ -7II 2 



21 ir^ >r^ {\i^ - A,, J Xj^Xj^^j 



1 2X1 _ '^32 



-7OO 

M 

'^32' 



-Til 1 °g" ^10 

>-'21 + ,^ >^22 



-7OO 



^11, 26;;^AiA2A3 .,-^ (Ai^ - Aij2(^^ _^2)2AiiAj, 

i721 ' / 2 OjiOi2''l''2 

^i''^ j^i2,i3,23 n4e/nj=i(A« + Mj) 



(4.43) 



and 



27166""^"^ 



ao oi + a2 + 03 
2m e""^" = '^°° '^^^ '^"^° oi + 02 + 03 



-Tio -Tio 7-10 3 2, ' 

>^iWoo >-'ii v^v^ Ai 

2=1 J = l "* 

-Til -Tio 

o„,„-2;4 _ ^10^^11 

^"-4^ — j-oi j-oi ' (4 44) 



2?TT,3e' 



-Til -tOI 

-X3 _ '^11 '^21 



10 -TlO' 



2n2e--^ - ^^2W2^ 



-TiU -Til 

>^22 >^11 

-Til -Tio 

■^21 '-'2: 
-tOI -tO 
>^21 >^32 



tOI -tOI ' 



2"^ie-^^ - y^ + 26:, - YYY^ + 2^° 

ij o'22 ^1-^2-^3 >-^22 



.1 ^ MJ2¥ ^ 26* ^ M1M2 >72Y 
LJ722 A1A2A3J22 

(The last few right-hand sides are too large to write in expanded form here, but 
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explicit expressions for the sums jZ^f„ are written out in Example A. 2.) 

Remark 4.12. Using the Lax pairs for the Geng-Xue equation, it is not difficult 
to show (details will be published elsewhere) that the peakon ODEs (1-13) induce 
the following time dependence for the spectral variables: 

A. = 0, d, = ^, /i,=0, bj^^, 6oo=0, b*^^0. (4.45) 

Ai flj 

This means that the formulas in Corollary 4.5 give the solution to the peakon 
ODEs (1.13) in the interlacing case, if we let the variables {\i, ^j,boo,bl^} be 
constant, and let {ai,bj} have the time dependence 

a,{t) = a,(0) e*/^% 6,(i) = 6,(0) e*/^^ (4.46) 

and the coefficients derived in Section B.3 are constants of motion. 

In particular, (4.41) and (4.42) give the solution to the 2 + 2 interlacing 
peakon ODEs 

xi = (mi + m3Ei3){n2Ei2 + n^Eii), 
±2 = {miEi2 + m2,E23){n2 + 714^24), 

is = (mi£;i3 + mg) (712^^23 + '^4£'34), 

±4, = (mi £'14 + m3^34)(n2-E24 + "4), 

TTZi 

= (mi + m3£;i3)(n2£i2 + ^4^14) - 2m3Ei3{n2Ei2 + n^Eu), 

mi (4.47) 

— = {~miEi2 + mj,E2i){n2 + niE2i) - 2{miEi2 + m^E23)n^E2A, 
n2 

Tfl'x 

= (miE'ia + m-i){-n2E2-i + ^4^:34) + 2miEi3{n2E23 + ?^4£'34), 

ma 

— = (-mi£;i4 - m3£'34)(n2i?24 + n^) + 2{miEu + m3£'34)n2i?24, 
71.4 

where Eij = g-l^i^^jl = g^i-xj jt^j. j ^ j^ ^^^ (4.43) and (4.44) give the solution 
to the corresponding ODEs for the case K ^ 3. Likewise, (4.55) below gives the 
solution to the 1 + 1 peakon ODEs 

mi 712 

±1 = 2:2 = — = = min2£'i2. (4.48) 

mi 712 

However, in this last case the equations are rather trivial, and all the heavy 
machinery is not really required. Indeed, m.in2-Ei2 is a constant of motion, so 
direct integration gives xi{t) — xi(0) + ct, X2{t) = 0:2(0) + ci, mi(i) = mi(0)e'^*, 
7i2(i) =7i2(0)e-=*, where c = mi(0) 712(0) e="i(°)-"=^(°). 

4.3 The case K = 1 

As already mentioned in Remark 3.2, the case X = 1 is degenerate. We have 



5(A) - L2(A) 



Li(A) 



io(A) 



1 - Xgihilo hi 

-Xgilo 1 

^A^5i/ii^o^2 — 2A —\h1l2 
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and 



Six) = L^iX) 



LoiX) 



1 51 

-Xhi{ln + li) 1-Xgihili hi \ , 
-2A -A.gi(;i + /2) 1 



and it follows (cf. (3.40)) that 

^21 (A) 



where 



and 



where 



W{X) = - 



Ai 



51^0 



ai 



SsiiX) Xgihilol2-2 A - Ai ' 



gihikh 



ai 



hih' 



^31 (A) 2 

hiik + h) hi{2~l2) 



Moreover (cf. (3.52)), 



where 



We also have 



531(A) 2 

,91(^1 + ^2) .91(2-/0) 



W*{X) 



532(A) 



hih 



S^iiX) Xgihilol2~2 A-Ai' 

where 

1 

oi = —r, 
gih 

so that (cf. Theorem 3.15) 

1 1 Ai 

aiQi = -; — — • — = -— . 

iT-ih gih 2 
From these equations it follows that 

r = aiOoo+ o' 

1 , 1 

J- = aiboc. + -, 



51 = Ko + 



hi = 60 



1 
2^1 

1 
2ai 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



(4.53) 



(4.54) 
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h^^' 


= fll&oo, 


lp-2xi 


* 7 * 


27126-^2 


1 

ai 
1 



Mapping back to the real line using — 1 + /q = yi = tanh 2^1 , 1 — ^2 = 2/2 = tanh 0:2 , 
gi — 2mi coshii and hi — 2n2 cosh 2:2, we get 



(4.55) 



277116"^ 

which recovers Xi, X2, rrii, n2 from the spectral data Ai, ai, boo, b'^ (and 
a\ = (2ai)-iAi). 

It is clear that all the spectral variables are positive if mi and n2 are positive, 
but there is an additional constraint (not present for K > 2): the ordering 
requirement xi < X2 is fulfilled if and only if 

1 < e2(-^--i) = Aaia\boobl, = 2\iboobl,. (4.56) 

In the terminology of Definition 4.7, the range of the forward spectral map S for 
K = \ is not all of TZ, but only the subset where A 1600^^ > \- 

5 Concluding remarks 

In this paper we have studied a third order non-selfadjoint boundary value 
problem coming from the Lax pair(s) of the nonlinear integrable PDE (1.2) 
put forward by Geng and Xue [14]. We have given a complete solution of the 
forward and inverse spectral problems in the case of two positive interlacing 
discrete measures. The main motivation for this is the explicit construction of 
peakons, a special class of weak solutions to the PDE; more details about this 
will be given in a separate paper. This inverse problem is closely related to 
the inverse problems for the discrete cubic string appearing in connection with 
peakon solutions to the Degasperis-Procesi equation (1.7), and for the discrete 
dual cubic string playing the corresponding role for Novikov's equation (1.5) (see 
[23] and [17], respectively), but it has the interesting new feature of involving 
two Lax pairs and two independent spectral measures. 



A Cauchy biorthogonal polynomials 

The theory of Cauchy biorthogonal polynomials, developed by Bertola, Gekhtman 
and Szmigielski [6, 5, 4, 7], provides a conceptual framework for understanding the 
approximation problems and determinants that appear in this paper. In Sections 
A.l and A. 2 below, we recall a few of the basic definitions and properties, just to 
give a flavour of the theory and put our results in a wider context. Section A. 3 is 
the crucial one for the purpose of this paper; it contains determinant evaluations 
(and also defines notation) used in the main text. 
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A.l Definitions 

Let a and /3 be measures on the positive real axis, with finite moments 

ak^ Jx'^daix), /3fc = y"/d/3(y), (A.l) 

and finite bimoments with respect to the Cauchy kernel l/(x + y), 



x + y 



(A.2) 



According to (A. 17) below, the matrix (labYal^Lo ^^^ positive determinant -D„ 
for every n, provided that a and /3 have infinitely many points of support. Then 
there are unique polynomials (Pn(a;))J^Q and (g„(y))J^o ^^*^^ that 

(i) degp„ = deg g„ = n for all n, 
(ii) the biorthogonality condition 



da{x)dj3{y) = 5ij 
X -\- y 

holds for all i and j (where 5ij is the Kronecker delta). 



(A.3) 



(iii) for each n, the leading coefficient of p„ is positive and equal to the leading 
coefficient of q„ . 

These polynomials are given by the determinantal formulas 



Pn{x) 



Qn{y) = 





^00 


hi ■ 


h,n-l 


1 


1 


ho 

I71-1.0 


hi ■ 
In-1.1 ■ 


Il,n~l 
In-l,n-l 


X 


\/DnDn+l 


a;"- 




InO 


hi ■ 


h,n-l 


x" 




Iqo 


hi •• 


h,n-l 


lon 


1 


ho 

In-1,0 


hi •• 
hi-1.1 ■ ■ 


Il,n-1 
hi-l.n-1 


hn 


VDnDn+1 


h-1 




1 


y 


yn~l 


yn 



(A.4) 



(A.5) 



If either a or /? (or both) is a discrete measure, the determinant _D„ will be zero 
for all sufficiently large n, and then there will only be finitely many biorthogonal 
polynomials; cf. Section A.4. 



A.2 Four-term recursion 

One basic property of the Cauchy kernel which underlies much of the theory is 
the following: 



h+l,b + Ia,b+1 — 



^a+lyb 



da{x)d/3{y) 



x + y 
x''y^da{x)dl3{y) 



j.ayb+1 



da{x)dj3{y) 



x + y 
x^-daix) f y'd^iy) = aaPb- (A.6) 
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For example, ii X ~ {Xij)i,j>o and Y = {Yij)i,j>o are the semi-infinite Hes- 
senberg matrices (lower triangular plus an extra diagonal above the main one) 
defined by 



X I Pi(^) 



X 



^Po(x)^ 
Pi{x) 



^90(2/)^ 



^go(y)^ 



91(2/) \ =Y \ iT^iy) 



then it is straightforward to show that (A. 6) implies 



Xn 



■Y„ 



TrnVn 



where the numbers 



Pn{x)da{x) 



^/D^ 



n+l 



^00 



nO 



hi 
hi 

hi 






(A.7) 



(A., 



are positive by (A. 19) below, and similarly for rjm = J qm{y)dp{y). Since the 
matrix 



f-^o 







L:= 



V 







-1 

3 


...\ 




■ ■•/ 



(A.9) 



kills the vector tt = (ttq, tti, . . . )^, we have = L-rrt]'^ = L{X+Y^) = LX+LY'^. 
If M[a_b] denotes the set of matrices which are zero outside of the band of diagonals 
number a to 6 inclusive (with the main diagonal as number zero, and subdiagonals 
labelled by negative numbers), then LX G -^[0,1] ' Af[-oo,i] = -^[-00,2] and 
LY^ e M[n.ii • Aff_i „oi = -A-^f-i.ooi- But since their sum is zero, it follows that 



[-1,00] ■ 



LX and LY are both in M[_i 2] (four-banded). Thus 



^Po{x)^ 
xL I Pi (2^) 



LX 



fpoix)^ 
Pi{x) 

V : 



(A.IO) 



with L G -M"[o_i] and LX G M[_i 2], is a four-term recurrence satisfied by the 
polynomials Pn- 

The same argument applied to X'^ + Y ~ rjiv^ shows that the polynomials g„ 
satisfy a corresponding four-term recurrence 



yL 1 9i(y) =Lr 9i(y) 



(A.ll) 



where L is like L except for n being replaced by rj. 

Further generalizations of familiar properties from the theory of classical 
orthogonal polynomials include Christoffel-Darboux-like identities, interlacing 
of zeros, characterization by Hermite-Pade and Riemann-Hilbert problems, and 
connections to random matrix models; see [6, 5, 4, 7] for more information. 
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A. 3 Determinant identities 

For X — {xi, . . . , Xn) G R", let 

A(a;) = A(a;i , . . . , a;„) = ]^(a;j - x^) 



i<j 



and 



T{x) = r(a;i,.. .,Xn) = Y\_(^i + ^j)^ 



(A. 12) 



(A.13) 



*<j 



where the right-hand sides are interpreted as 1 (empty products) if n = or 
71 = 1. Moreover, for x £ R" and y E R™, let 



4=1 J=l 



(A. 14) 



Finally, let (T„ be the sector in R" defined by the inequalities < Xi < • • • < a;„. 
With this notation in place, we define 

J:L = / ^, ^^^ ^-da-{x)dr{y), (A.15) 



for n and m positive. We also consider the degenerate cases 



J:^= A{xf(Y[xXda-{x) (n>0), 

JoZ^ [ A{yr(flyXdr{y) (m>0), ^^'^'^^ 

>-'oo — ^■ 

Note that J{q = a^ and J^l — /3s . 

When a and /3 are discrete measures, the integrals J^^^ reduce to sums; see 
Section A. 4 below. 

Many types of determinants involving the bimoments lab — J J i;;^da{x)dl3{y) 
can be evaluated in terms of such integrals; these formulas are similar in spirit to 
Heine's formula for Hankel determinants of moments ak = J x da{x), encoun- 
tered in the theory of orthogonal polynomials: 



00 



det(a,+,)^-io= / A{xrda-{x) = X^> 

Here we have collected a few such formulas (all formulated for determinants of 
size n X n). To begin with, specializing Theorem 2.1 in [6] to the case of the 
Cauchy kernel K{x,y) = l/{x + y), we get the most basic bimoment determinant 
identity, 

-^00 • ■ ■ lQ,n~l 



D„ = det(/,, ^"-1 



^ijlij=0 



In- 



1,0 



In- 






00 
nn' 



(A.17) 
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Applying (A. 17) with the measure x^ da(x) in place of da{x) and y'' dp{y) in 
place of d(3{y) gives 

det(/,+,,,+,)^7io = J--. (A.18) 

Proposition 3.1 in [6] says that 

^00 • ■ ■ Io,n-2 



In-i,a 



ao 



In-l.n-2 O^n-l 



*-^ n,n—l • 



(A.19) 



(For n — I, the left-hand side should be read as the 1x1 determinant with the 
single entry ao; this agrees with J^iq = ao-) By the same trick, we find from 
this that 



-'r+n— l.s • ■ • -'r+n— l,s+n— 2 ^r-\-7i—l 

We also need the following identity: 

do ho ■ ■ ■ Il,n-2 



J-rs 



(A.20) 



*-^ n,n—l 



(A.21) 



Ckn-l InO ■ ■ ■ In,n-2 

Proof of (A.21). From (A.20) we have 

-^01 • ■ ■ lo.n-l 



rrOl 



ao 



In~l.l 



-'n— l.n— 1 Ctn— 1 



We rewrite the bimoments as Ijk = aj/3j,_i — /j+i^fc-i (using (A. 6)), and then 
subtract /3k~i times the last column from the other columns k — 1, . . . ,n — 1. 
This transforms the determinant into 



'10 



-InO 



-Il,n-2 Q^O 



-In,n-2 "n-l 



without changing its value. Now move the column of a's to the left; on its way, 
it passes each of the other columns, thereby cancelling all the minus signs. D 

Other useful formulas follow from the Desnanot-Jacobi identity, also known 
as Lewis Carroll's identity (or as a special case of Sylvester's identity [12, Section 
II. 3]): if X is an 71 X n determinant (with n > 2), then 



v\^ vn7i "v^ll V^^ Y^^^^ 



(A.22) 



where X^^ is the subdeterminant of X obtained by removing row i and column j, 
and where Y — (X"")^^ is the "central" subdeterminant of X obtained by 
removing the first and last row as well as the first and last column (for n = 2, we 
take y = 1 by definition). For example, applying this identity to the bimoment 
determinant (A.18) (of size n + 1 instead of n) gives 



^rs ,T-r+l,s+l ^rs q 



7's //r+ljS+l 
nn 



^ nn ^ nn ' 



(A.23) 
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and from (A.21) we get 

(for n> 1, in both cases). 

Another identity, which holds for arbitrary z^, Wi, X^j, is 



(A.24) 



-S-n ^ n 1 



W2 -^^21 



-^l,n-l 
X2,n-1 

^n,n— 1 



W2 



-^^21 



-'^nl 



l,n-l 
-'^2.n-l 



x„ 



Wn Xnl 

Wn+1 Xn+IA 

Xn 



Zl 
Z2 



'-n+l 



Xn 

^21 



Xnl 



Xl^n-1 
X2,n-1 

^n,n— 1 



-'^2n 



2^1 
22 



-''^21 

u;2 



Xl.n-l 
X2,n-1 

^n,n— 1 
^n+l,ri-l 

Xi^n-1 
X2,n~l 



X. 



n.n— 1 



X 



n+l,n— 1 



Xln 
X2n 

Xn+l,n+l 
X2n 

n+l,n+l 



X. 



Xn 



2n+l Wn+1 



Xnl 
Xn+1,1 



X2,n-1 
-^n.n—1 



(A.25) 



Indeed, the coefficients of 2„+i on the right-hand side cancel, and the coefficients 
of the other variables Zj on both sides agree, which can be seen by applying the 
Desnanot-Jacobi identity to the second determinant on the left-hand side with 
its jth row moved to the top. 

In the text - see equation (4.5) - we encounter the n x n determinant 



ICn 



foo + 2 

ho 



ho 
ho 

ho 



h,n-2 
h,n-2 

I n.n— 2 



(A.26) 



which satisfies the recurrence 

ICn+l ICn 



-tui -tui 



-710/ -TOO I 1 -Til \ 

>^rmWn+l,n ^ 2'-^n,n—l) 



(A.27) 



Proof of (A.27). By taking Zi = ai-i, wi = Iio + \5iQ and Xij — h.j-i in (A.25), 
and using (A.18) and (A.21), we find that J°]n-i^n+i ^ ICnJn\i,n " Jnl^, 
where 



Z = 



ao h 
ai 



00 -^ 2 
h 



'10 



ho 
ho 



hi 
hi 



l,n-2 
2,n-2 



InO 



h 



^n+1,0 -^n+1,1 

Using again the trick of rewriting the bimoments (except in column 2) as 
Ij+i.k = Oijfik — Ij,k+ii subtracting Pk-s times the first column from column k 
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(for k — 3, . . . ,n), and moving the first column to tlic right, we see that 



Z = - 



-^10 hi Il2 



Iq.ti-I 0!o 

Il.7l-1 Oil 



InO In,l In,2 ■ ■ ■ In,n-1 Q^ji 

Wri+l,ri ^^ 2'-'n,n—ll- 

(The last equality follows from (A. 20).) Consequently, 

-701 r — r -701 -I- ^Wf^oo I 1 711 ^ 

which is equivalent to (A. 27). n 

Remark A.l. It was shown in [23, Lemma 4.10] that when a = /3, the factor- 
ization 

1 / /• A(a;)2 






T{x) 



da"(a;) 



(A.28) 



holds. We are not aware of anything similar in the general case with a ^ (3. 

A. 4 The discrete case 

Consider next the bimoments lab and the integrals J^m defined by (A. 15) in the 
case when a and j3 are discrete measures, say 



i=i 3=1 



(A.29) 



(In the setup in the main text, we have A ~ K and B = K — 1.) Then the 
bimoments become 

lab - // ^da{x)dP{v) = J2J2 Y^ "'^^■' (^-30) 

J J x + y i=i j=i '^^ ^ ^^3 



and likewise the integrals JTJ^,^^ turn into sums: 

^6C^')-^e(lfl) 



(A.31) 



Here ('^') denotes the set of n-element subsets / = {zi < ^2 < • • • < in\ of the 
integer interval [A] = {1, 2, ... , A\, and similarly for ( „j). Morever, 



and 



is/ jGJ 



VI.,-- ^?^'- 



'/J 



/J 



(A.32) 



(A.33) 
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where we use the shorthand notation 

A2 = A(A,,,...,A,J2= l[{Xa-Xb? 



a.bel 
a<b 



A^J^A{^lJ,,...,^i,J^^ lli^ia-f^b)' 



a,be.J 
a<b 



(A.34) 



Tij —T{\i^,.. . ,Xi^;fij^,.. .,fij^) — 11 (Xi + fij). 

ieijeJ 

For later use, we also introduce the symbol 

A?,/.- n (^n-A.J^ (A.35) 

and similarly for Aj^j^- Empty products (as in A^ when / is a singleton or the 
empty set) are taken to be 1 by definition. When needed for the sake of clarity, 
we will write ^/,,/ instead of ^/j, etc. 

For (positive) measures on the positive real line (a^, A^, bj, iij positive), we 
thus have J^^^ > for < n < A and < m < B, otherwise J7^^ = 0. 

Example A. 2. Below we have listed the nonzero jT^^j in the case A = 3 and 
B — 2 (the more general sum J7^^ is obtained by replacing each a^ and bj in 
Jnm by A[aj and ^i^bj, respectively): 

Jio = ai + 02 + 03, 

Jw = (-^1 - -^2)^0102 + (Ai - \3)^aia3 + (A2 - \3fa2a3, 



J^^ = (Ai - A2)^(Ai - A3)^(A2 - A3)^aia2a3, 

-7UO _ T 



Jo°i° = bi + b2, 



''01 

rOO 



T — ; aih + T — ; 0261 + T ; 0361 

Al + /il A2 + Ail A3 + /ii 

-aib2+- ; 0262 + T ; 0362, 



Al + fl2 X2 + H2 A3 + fl2 

^00 (Al - A2) (Al - A2) 

'^si = TV": ^TV": raia2foi + 77—- tt^-— 7010262 

(Al +^i)(A2 + Afi) (Al +M2)(A2 + M2) 

(Al - Xs? , , (Ai-A3)^ , 

[Xi+ ^l){X3 + fli) (Al +^2)(A3 + M2) 

(A2 - A3) (A2 - A3) 

+ TT , TTT ^ 7020301 + TT , TTT , 7020362, 

(A2 +Ail)(A2 +Ml) (A2 +A*2)(A2 + M2) 

J3^i> - (^r'^^yr'^ll'^-'^^iwi 



2 



(Al + Aii)(A2 + Mi)(A3 + Ml) 
(Ai-A2)^(Ai-A3)^(A2-A3)^ 

(Al + /i2)(A2 + A*2)(A3 + H2) 



54 



01020362, 



^00 (a^i - Ma)^ , , , (mi - M2)^ , , 

(Xi + fli)(Xi + fl2) (A2+Ail)(A2+M2) 



-03616; 



'2, 



{X3 + Hi){X3+ fl2) 

TOO _ (Ai - X2)Hfii - M2)' , , 

(Ai + Ml)(A2 + A^l)(Al + /i2)(A2 + ^i2) 
(Ai - A3)2(a.1 - /i2)' 



(Ai + Aii)(A3 + A*i)(Ai + At2)(A3 + H2) 

{X2-x^f{^l^-^l2f 



aia3bib2 
02036162, 



(A2 + /il)(A3 + fIi){X2 + M2)(A3 + fJ-2) 

joo ^ (Ai - A2)^(Ai - A3)^(A2 - A3)^(a.i - (^2)' q^a^g^b^b., 

^^ (Ai + /Xi)(A2 +Ail)(A3 + /il)(Ai +^2)(A2 + Ai2)(A3 +A*2) 

Lemma A. 3. Let (i/*)^^*,! denote the integral iZ^,^„ evaluated using the measures 

A B 

a*=^a*(5A, and /3*==^6*5^^. (A.36) 

i=i j=i 

in place of a and j3, and suppose that the Jour measures are related as in 
Theorem 3.15: 






Then 

2?n— n+s— 1 



^ A~n,B-m 



^- ' \^=i / '^^^ (A.38) 



'-' A-n.B-m 



On+7n rr^n—2n+l—r,n—2m+l — s 
^ J A 



' AB 



Proof. This is a fairly straightforward computation. To begin with, let L 
Y\i=i Aj and M — ni=i Mi; a-^id write (A. 37) as 



OfeOfc 





L^T{k},iB] 




6, 


,* ^'r[^],^,j 




A 


' 


cOfc B 


2A* 
;her 


.M Y[{X, - 
efore have 


A.)^ 






n* - 




1 / 

a7U 


2M ) 


^/.[B] 


Oj - 
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(The factor in front is A/ = Yliei ^i-) ^^ ^^'^ denominator, the factor (Ap — A^)^ 
will appear twice if p and q both belong to /, once if one of them does, and not 
at all if both belong to [A] \ I. Thus we can write 



fl/flj = 


1 
^ a7 




1 ^'i^u ■ 


bjb:, = 


1 


\2Lj 


' A2 A^ 



Similarly, 



Putting this into the definition of {J*Ynm gives 



Now note that 



which implies 



A|A5 ^_ / L^y IY[BlA|<jy 



EE ie^; 



Yij "' Xiai \2M J A2A2 



2 

X 

[A] 



.2 



, 1 (M^Y^[A],.J^^B]\.J 
M.7 



^^JhJ\2LJ A^.Af^, 



r/,[B] — r/,jr/,[B]\j, 
^[A],.j = r/,jr[A]\/,.7, 
r[A],[B] =- ^ i,.j^[A]\i.j^ i.\B]\j^[A]\im\j^ 



^I,[B] ^[AIJ _ ^[AUB] 



^I,J ^IA]\L[B]\J 
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Thus, 



W )nm 7 J 7 J 



I J ^ ^ ■' 

1 [A],[B] ^[A]\/^[S]\J 

X 



X 






V J vU[A]\/y U[i3]\J 2"+ 

Q[A]\jfe[B]\J ^[A]\J,[B]' 

J 

T 2n—m-\-r—l A/r2m— n+s— 1 

■X 



H XI *[^]U,[S]\./ i\A]\iy^'' a[A]\I it^lB]\j)^^' b[B]\J 
I J 

T 2n~ni+r—l )\/r2T??.— n+s— 1 

X"^ V^ T \ 1 — r 1 — si 

= 2"+" 700 ^ ^ "^uv^u ^uf^v ov 

T 2n—m+r—l i\/f2in—n+s—l 
_ ^ f*^' rrl-r,l-s 

~ On+m -TOO '-'A~n,B~'m^ 

as claimed. D 

The following lemma is needed in the proof of Lemma A. 4. 

Lemma A. 4. Suppose (as in the main text) that the number of point masses in 
a and /? are K and K — 1: 

K K-1 

Then the quantities J'^m satisfy 

Jjj-iJj-ij-i ~ Jjj Jj-ij-2 > Uj J — 2, . . . , A! — 1, (A.d9) 

and 

J™J^:i_i-J^i,J^i,_i>0, j^l,...,K~l. (A.40) 

Proof. In (A. 39) we let j = m+ I for convenience {1 < m < K — 2), and expand 
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the left-hand side using (A. 31): 

-tOO n-ll _ /xOO -xll 

^ m+l ^m*-^ m7n *-^r??.+ l.rn+l*-^m,m — 1 

= I XI Yl '^AC «A ^c 1 I XI X *SD As Ai£, ag b^ j 

- I X X ^-40 «A ^c I I X X *SD As Md Ob bo j 

= X X X X *Ac*BD -^s Md flA as ^C ^D 

AG^^'i) bg(I^I) ce(i^^ii) ceC^-^O 

~ H X X X ^AC^BD^sMDaAas^c^D- 

MiVi) ^e('^0 ceC^;;0 ce(i^:;i) 

(A.41) 

Let us denote the suniniand by 

f{A, B, C, D) = '^Ac'^BD ^B Md a a ag be b^ 

for siniphcity. 

Choosing two subsets A and B of a set is equivalent to first choosing R — AOB 
and then choosing two disjoint sets X — A\{Ar)B) and Y = B\{Ar)B) among 
the remaining elements. If | A| = m + 1 and \B\ — m, then 

\AnB\ = m-k, \A\{AnB)\ = k + l, \B\{AnB)\ = k, 

for some k e {0, 1, . . . , m}. Thus 

m 

X E fiA,B,C,D)^Yl E E fiR + X,R + Y,C,D), 

xny=0 

where we write R + X rather than RU X, in order to indicate that it is a union 
of disjoint sets. (If K — (m — fc) < 2fc + 1, then the innermost sum is empty.) 
With similar rewriting for C and D, (A.41) becomes 

rn m~l 

EE E E E E f{R+x,R+Y,s+z,s+w) 

ye(m\«) wg(1^71\^) 

m m— 1 

-EE E E E E /(i?+^,i?+>^,^+^,^+M^). 

fe=o /=o ^g( i.Ky gg(iK-i,)^g([^^\«) ze(i--i\-) 

re(i-ix«)H'6(iVl^-) 
jfnF=0 znVK=0 

(A.42) 
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(The innermost sum on the second hne is empty when I = 0.) Now, since 

a2 X2 a2 X2 

^ R+X.S+Z^ R+Y.S+W — -p p 

-■- R+X,S+Z -•- R+Y,S+W 

A 2 A 2 A 2 X2X2X2 A 2 A 2 A 2 a:2a.2 a:2 
^ ^R^X^RX^S^Z^SZ ^R^Y^RY^S^W^SW 

^RS^RZ^XS^XZ ^RS^RW^YS^YW 

A 4 A 2 A^A^ 

— AAj^AAyAA^^VF^ XW'^ YZ -r^ F ' 

-•^ fl+X+F,S+Z+W^ RS 

we have 

/(i? + X, i? + r, 5 + z, s* + VF) 

= ^ R+X,S+z'^ R+Y,S+W ^R+Y f^S+W ^^R+X '^R+Y ^S+Z ^S+W 
= {^X^Y^l^W^XW^YZ^Yf^w)^ 

/a^a^ a^a^ \ 

/ ^R^R,X+Y^S^S,Z+W , 2 ;,2 , \ 

~^ r '^RfJ'S '^R'^X+Y "S "Z+W y 

\ '^ R+X+Y,S+Z+W'^ RS J 

where the first factor depends on the sets X,Y,Z and W individually, while the 
second factor involves only their unions U = X + Y and V — Z + W. Therefore 
we can write (A. 42) as 

fl.e(,rj,) 5e(if„-i) r/e(S)f) ye('--l^-^ V ^R+u,s+v^RS 
•^flMsa|ay&sV Yl Yl ^x^Y^l^w'^xw'^YZ>' 

\X+Y=U Z+W=V 




Y^^W 



|X|=fe+l \Z\=l 
\Y\=k \W\=l 



(A.43) 



^2 a2 a2a2 

X+Y=U Z+W=V 

|X|=fc+l \Z\=l+l 

\Y\=k \W\=l-l 



A^Ay^Z^W^XW^YZ'^Yt^W 



This is in fact positive (which is what we wanted to prove), because of the 
following identity for the expression in brackets: if \U\ = 2fc + 1 and \V\ = 21, 
then 



7 J / , ^X^Y^Z^vi^ XW^YZ^Yi^W 



X^rY^U Z+W=V 
\X\=k+l \Z\=l 
\Y\=k \W\=l 



/ , 7 , ^X^Y^Z^W^ xw^yz^yI^w 

x+Y=u z+w=v (A. 44) 

\X\=k + l \Z\=l + l 
\Y\=k \W\=l-l 






7 , '-^X'-^Y'-^Z'-^W^XZ^YW^Y^^W■ 
X+Y=U Z+W=V 
\X\=k+l \Z\=l 
\Y\=k \W\=l 
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(Note the change from ^ xw^yz °^ ^^^ ^^^^ *o ^ xz^yw '^^ ^^^ right.) 

To prove (A. 44), we can take U — \2k-\-l\ and V = [2/] without loss of general- 
ity. If ^ = 0, the identity is trivial, since both sides reduce to '^x+y=u ^x^y-^f- 
The rest of the proof concerns the case Z > 1. 

When fc = and I — 1, both sides reduce to Ai(/ii + /X2) + 2fii^2- For 
fixed {k,l) with fc > 1, evaluation at X2k = X2k+i — c gives on both sides 
2c J^j^j^ (Xi — c)^ ni=i('^ + Mj) times the corresponding expression with fc — 1 
instead of fc. Provided that the identity for (fc — 1,Z) is true, our {k,l) identity 
therefore holds when A2fc = X2k+i, and in fact (because of the symmetry) 
whenever any two Xi are equal. This implies that the difference between the 
left-hand side and the right-hand side is divisible by A^. (Any polynomial p in 
the variables A^ which vanishes whenever two A^ are equal is divisible by A[/. 
If in addition p is a symmetric polynomial, then p/Au is antisymmetric and 
therefore again vanishes whenever two A^ are equal; hence p/Au is divisible by 
A[/.) Considered as polynomials in Ai, the difference has degree 2fc + I and Afj 
has degree 4fc. If I < 2fc, then 2k + I < 4fc, and in this case we can conclude that 
the difference must be identically zero. To summarize: if the (fc — 1, Z) identity 
is true and I < 2fc, then the (fc, /) identity is also true. 

Similarly, for fixed (fc, /) with / > 2, evaluation at fi2i-i = M2i = c gives 
2c Jlj^-^ (Ai + c) Y[j=i (Mj ~ c)^ times the corresponding identity with I — 1 
instead of L As polynomials in /ii, the difference between the left-hand side 
and the right-hand side has degree 21 + k and Ay has degree Al — 2. The same 
argument as above shows that if the (fc, I — 1) identity is true and k <2l ~ 2, 
then the (fc, I) identity is also true. 

Since any integer pair (fc,/) with fc > and I > 1, except (fc, /) — (0,1), 
satisfies at least one of the inequalities / < 2fc or fc < 2Z — 2, we can work our way 
down to the already proved base case (0, 1) from any other (fc,/) by decreasing 
either fc or Z by one in each step. This concludes the proof of (A. 44), and thereby 
(A. 39) is also proved. 

The proof of (A. 40) is similar: the left-hand side expands to 



EE J2 J2 J2 T. fiR+x,R+Y,s+z,s+w) 

/c=o 1=1 ^g( M) se(>--i) x^i^-i}-) zg(1-7)\-) 
ye([-i\«) we(i--l\-) 

j j 

'EE E E E E f{R+x,R+Y,s+z,s+w), 

fc=o 1=1 «g( [-1 ) sg([^_-ii) X6(i-|\«) ze(i--i\-) 
ye(i-l\f ) we(i--]\-) 
xny=0 znVK=0 

(A.45) 
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which equals 



EE E E E E r % X 

\X+Y=U Z+W=V 

\X\=k \Z\=l 
\Y\=k \W\=l-l 



(A.46) 



E E ^^^y^l^'v^r^M/ry^A^M,^) J. 



X+Y 

\X\=k+l \Z\=l 

|y|=fe-i \w\=i-i 
which is positive, since for \U\ ~2k and |V^| = 2/ — 1 the identity 

/ , / , ^x^y^z^wrxn'-TyzAy/iT^ 
x+y=;7z+w=y 

\Y\=k \W\=l-l 

— 2^ 2^ ^X^Y^Z^W^XW^YZ'^Yt^W 

x+Y=u z+w=v (A. 47) 

\X\=k+l \Z\=l 

|y|=fc-i \w\=i~i 

= 2_^ Z^ ^X^Y^Z^W^XZ^YW^Yi^W 

X+Y=U Z+W=V 
\X\=k \Z\=l 
\Y\=k \W\=l-l 

holds; it is proved using the same technique as above. n 

B The forward spectral problem on the real line 

Consider the Lax equations (1.3a) and (1.4a) in the interlacing discrete case (3.1). 
In what follows, we will start from scratch and analyze these equations directly 
on real line, without passing to the finite interval (—1, 1) via the transformation 
(2.1). As we will see, this leads in a natural way to the definition of certain 
polynomials {^fe(A), -Bfe(A), Cfe(A)}j^g. Since the two approaches are equivalent 
(cf. Remark 2.1), these polynomials are of course related to the quantitites defined 
in Section 3. Before delving into the details, let us just state these relations, for 
the sake of comparison. 

In the interval yk < y < yu+i, the wave function ^(j/; A) is given by 

^0i(2/;A)\ / Ak{\)-\Ck{\) \ 

02(2/; A) = -2ABfe(A) . (B.l) 

Mv^ A)/ V-A(i + y)M>) + A2(i - y)Ck{X)l 

Hence, letting y — > 1^ we obtain (with {A, B, C) as synonyms for {An, Bn, Cn)) 

Vi(l;A)\ /A(A)-AC(A)\ 

(^2(1; A) = -2Ai?(A) . (B.2) 

,03(1; A)/ V -2A^(A) J 
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As for the spectra and the Weyl functions, we have 

K , , X K-l 



A(A).n(i-^)^ 1(A). n(i- A), (B.3) 

'^^'^-"0^(iyA)--A(A)' ^(')--0^(M)-2A-A(A)' (^-^^ 

l^(A) = -&iilA = _|W, Z(A) = -tliil^^i--2A). (B.5) 
03(1; A) A{\) 03(1; A) 2A A(A) 

The residues Oj, 6j and &oo are defined from the Weyl functions as before; see 
Theorem 3.10. One can also define similar polynomials corresponding to the 
adjoint Weyl functions, in order to define h*^ as in (3.52), but it is perhaps 
more convenient to define b*^ using the relations (3.56) and (3.58), which in this 
setting take the form 



bocblo = ( 



i=l 



K-1 



^ (fj (l--^2»-l,2») \ /t-T ^J-J 

2U^ A„:i?2Vi,2, Ay (l-i?|,-2,+i) 



K>2, 



K=l, 



(B.6) 



l2Ai£;2^' 
where Eij ~ e^'^'^^j' — e^'^^^ for i < j. 

B.l Setup 

To begin with, recall equation (1.3a), which determines ^(x; z): 

dxipi{x; z) = zn{x)'ip2{x; z) + ipsix; z), 

dxip2{x; z) = zm{x)i)z{x; z), (B.7) 

dxip3{x;z) = 'il;i{x;z), 

Away from the points Xk where the distributions m and n are supported, this 
reduces to 

dxipi^i's, dxip2^0, dxi'z^i'i, 

so 'ip2ix; z) is piecewise constant, and i^i{x; z) and ^'3(2;; '^) ars piecewise linear 
combinations of e^ and e^^. It is convenient to write this as 

'Mx;z)\ /Ake- + z^Cke--\ 

i{j2{x; z) \ = \ 2zBk , Xk<x<Xk+i, (B.8) 

^^3(0:^; z)J \Ake^ - z^Cke-^j 

where the coefficients {A^, Bk, C'^j^g may depend on z but not on x. (Here we 
set xo — —00 and xn+i — +00, so that the x axis splits into A^ + 1 intervals 
Xk < X < Xk+i numbered by fc == 0, 1, . . . , N .) Then the conditions (2.4) and 
(2.5) for ^(a;;z) at ±cx3 translate into 

Bo = Co = = An. Ao - 1, (B.9) 

respectively. So we impose (Ao,i?o,C'o) = (1,0,0) (i.e., ^(x; z) = (6^,0,6^^)-^ for 
X < xi) and investigate for which z the condition An{z) = is satisfied; the 
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corresponding values X — —z'^ will be the eigenvalues considered in the main 
text (cf. Remark 2.2). 

The pieces (B.8) are stitched together by evaluating equations (B.7) at 
the sites x — Xk- Since the Dirac delta is the distributional derivative of the 
Heaviside step function, a jump in ipi at Xk will give rise to a Dirac delta 
term 5x^ in dxtpi, whose coefficient must match that of the corresponding Dirac 
delta coming from to or n on the right-hand side of (B.7). Denoting jumps by 
[/(a;fc)] = fi^t) ~ /(^fc )i ^6 fi^^d ^t the odd-numbered sites x = Xk = X2a~-i 
(where m is supported) the jump conditions 

[ipi{xk;z)] = 0, 
[ip2{xk;z)] = 2zTOfcV3(a;fc), 
[■4;3{xk; z)] = 0, 

while at the even-numbered sites x = x^. ^ X2a (where n is supported) we get 

[ipi{xk;z)] =2znkip2{xk), 
[i!2{xk;z)\ = 0, 
[Tpz{xk;z)\ = 0. 

Upon expressing the left and right limits ipi{x^) using (B.8), these jump condi- 
tions translate into linear equations relating {Ak, Bk, Ck) to {Ak-i, -Bfc-i, Ck-i)- 
Solving for {Ak, Bk,Ck) yields 



'Ak\ ^ /Ak^i 

Bk I = Sk{—z ) I -B/c-i 

^CkJ \Ck-i, 



(B.IO) 



with the jump matrix Sk (not to be confused with the transition matrix S{X) 
defined by (2.12) and used in the main text) defined by 



Sk{X) = < 



' / 1 

nike^^ 1 Xruke^'^'' 
^00 1 




fc = 2a- 1, 



2a, 



(B.ll) 



for a = 1, . . . , A'. Starting with (^Oi ^o, C'o) — (li 0, 0) we obtain in the rightmost 
interval x > xn polynomials {Ajy, Bn, Cm) = {A{X),B{X),C{X)) in the variable 

1(A)) ^S2k{X)S2k-i{X)---S2{X)S^{X){}^. (B.12) 



C{\) 



Between the factors S2Ki^) and 5*1 (A) in the matrix product there are K — 1 
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S2a+lWS2aW = | m2a+ie"-+^ 1 

2n2ae^^» 

/O -2n2ae-^=^° 

A 2TO2a+i«2a(e"^"~"^"+^-e"^»+i-"^») m2a+ie-^^»+i I , (B.13) 
\0 

each such pair depending Hnearly on A. The factor S'i(A) (1,0,0)"^ does not 

depend on A, and therefore the vector (A(A), i?(A), C(A)) equals S2kW times 
a vector whose entries have degree iiT — 1 in A. Since A only appears in the top 
row of S'2i<-(A), we see that B{X) and C(A) are polynomials of degree K — I, 
while ^(A) is of degree K. We will name the coefficients in these polynomials as 
follows: 

A{X)^1-2X[A], + --- + {-2X)''[A]k, 

B{X) = [B]o - 2A[B]i + • • • + {-2X)''-^[B]k-i, (B.14) 

C(A) = [C]o - 2A[qi + • • • + (-2A)^-MC]k-i. 

These coefficients can be computed explicitly in terms of the positions Xk and 
the weights rrik and n^. by carefully studying what happens when multiplying 
out the matrix product S2k{^)S2K-iW ■ ■ ■ S2{X)Si{X) (f,0,0)-^. For example, 
using the abbreviation 

Eab = e^l"^""^"! (= e^"""^" when a < b) (B.f5) 

we have 

[A]i= ^ m^njE^j, 

l<i<j<N 

[A]k = min2Ei2 (1 - ^23) rn^ntEs^ (1 - £^45) m^riQE^^, ■ 

• • • (1 — -BAr_2.A'-l) ™Af--l'^Af£^JV-l,Af, 



(B.f6) 



[B]o = 5(0) - 



l<i<N 



(Recall that N — 2K. Note also that since m2a and n2a-i are zero, only the 
terms with i odd and j even contribute to the sums.) Later we will show a 
simpler way to read off all the coefficients in A{X); see (B.37) in Section B.3. 

For the second Lax equation (1.4a) things are similar, except that the roles 
of m and n are swapped. This leads to 

ip2{x;z)\^\ 2zBk |, Xk<x<Xk+i, (B.f7) 

2/ 



,^'3(3;; z)l \Ake'' - z^Cke 



and 




(B.18) 
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where 




fc = 2a- 1, 



SkW = { 



(B.19) 



/ 1 \ 

nfee^'= 1 XukC^^'' , k ~ 2a. 

IV 1 ; 

Starting again with {Aq, Bq,Co) — (1,0,0), we have in the rightmost interval 
iAN,BN,CN) = {A{X),B{X),C{X)), where 



f B(A) ) = S2k{\)S2K-iW ■ ■ ■ S2{X)Si{X) ( o) 



(B.20) 



Because of the asymmetry between m and n coming from the interlacing, the 
variable A appears in a shghtly different way here; in this case we have K pairs of 
factors S2aS2a~i, each of a similar form as the pairs S2a+iS2a that we computed 
earlier: 



S2a{\)S2a-l{\) = \ n2ae^^' 




-2m2a-ie-^^"-i 

A I 2n2am2a-i(e"^"-^-"^° - e"^"-"^-!) 7i2ae-"^° 





(B.21) 



From this we see that S2S1 (1,0,0)"^ is independent of A, so that the degrees 
of ^(A), B{X) and C(A) are at most K — 1. The leftmost pair S2KS2K-1 has 
no A in its bottom row, so C(A) is in fact only of degree K — 2. Naming the 
coefficients as 



(-2A)^'-MA] 



K-l, 



+ {-2X)''-'[B]k-i, 
+ (-2X)''-^[C]k-2, 



(B.22) 



A(A) = l-2A[A]i + - 
B{X) ^ [B]o - 2X[Bl 
C{X) = [C]o-2X[C]i- 
we have for example 

l<j<i<N 

[A]k-i = ^27713^23 (1 - Eji) nim^E^s (1 - E^q) nemrEej ■ 

• • • (1 — Ej^_^pf_2) nN^2'mN^iEN-2,N-l, 

[B], = B{0)= J2 ^J^"'^ 

l<j<N 

[B]k-i = [I]K-i«we""(l-S^_i,jv), 

where (as in (B.16)) only terms with i odd and j even contribute to the sums. 
See (B.40) in Section B.3 for an easy way to read off all the coefficients of A{X). 



(B.23) 
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B.2 Positivity and simplicity of the spectra 

By construction, the zeros of A{X) and ^(A) are exactly the nonzero eigenvalues 
Ai, . . . , Xk and fj,i, . . . , ^k-i treated in the main text, so the following theorem 
implies Theorem 3.8: 

Theorem B.l. // all nonzero weights m2a-i o-nd n2a o-^e positive, then the poly- 
nomials A{X) and A{X) have positive simple zeros Ai, . . . , Xk and /ii, . . . , Hk-i, 
respectively. 

Proof. We will rewrite the two spectral problems, (1.3a) with boundary con- 
ditions Bq = Cq = Q = Ajs! , and its twin (1.4a) with boundary conditions 
Bq = Cq — Q = Ajv, as matrix eigenvalue problems. (Recall that A and A are 
just aliases for Ajq and Ajv, respectively.) 

For the first problem, elimination of ■01 from (1.3a) gives dx'il^2 = zrmp^, and 
(9^ — l)'ip3 = znip2, which, considering the boundary conditions in the form (2.4), 
we can write as 

M^)^z My)dm{y), Mx) = -z -e-^^^-^^My) dn{y). (B.24) 



Evaluating the first equation at the even-numbered Xk and the second equation 
at the odd- numbered Xk, we find the 2K x 2K eigenvalue problem 



V'2,cvon\ ^ ( 



^' ■' "'V'")fc)- '--) 



where 

^2,cvcn = {'^2{X2),'4'2{X4), . . ■ ,^2{X2k)) , 
V'3,odd = (^'3(2^1), "03(2:3), • ■ • , i'3{x2K-l)) , 

I = K X K identity matrix, 

C ~ strictly lower triangular K x K matrix with Cij = 1 for i > j, 

£ = (e"l^^'-i^^^^l)5^i = (£;2^_i,2j)5=i (using the notation of (B.15)), 

M = diag(mi,m3,.. .,TO2k_i), 

JV = dia.g{n2,n4, . . . ,n2K)- 

Eliminating ^s.odd we can write this a,s a, K x K eigenvalue problem in terms of 
'02,ovcn alone: 

^2,cvon = 2A (/ + C)M£J^ ^2,cvcn (A - -z^). (B.26) 

As we've seen earlier, the eigenvalues are given precisely by the zeros of ^(A), 
and since ^4(0) = 1 we must therefore have 

A{X)^dct{I~2X{I + C)M£Af). (B.27) 

Now, for positive numbers {m2k-i,n2k}k=i, the matrix (/ + C)M.£M is oscilla- 
tory, since / -I- £ is nonsingular and totally nonnegative (being the path matrix 
for the planar network illustrated in Figure 5), and since M.£M is totally positive 
(£ being a submatrix of the totally positive matrix {Eij)f^^^). This implies that 
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>o KD 

>A >0 

>i .0 

-0 >Q 

Figure 5. A planar network (illustrated in the case K = 4) for which I + jC is 
the path matrix. What this means is that matrix entry (i, j) equals the number 
of paths from source node i on the left to sink node j on the right; in this case 
there is one path ii i > j and none if i < j. 



its eigenvalues, which up to an unimportant factor of 2 are the zeros of A, are 
positive and simple. (See, for example, our earlier papers [23, 17] for a summary 
of the relevant results from the theory of total positivity used here, and for 
further references.) 

For the second spectral problem we swap m and n and obtain 

/X /"OO -1 

My)dn{y), M^) = ~z -e-^--y^My)dm{y). (B.28) 

-OO J —OQ ^ 

These integral equations are evaluated the other way around (at odd-numbered 
and even- numbered Xk, respectively); this yields 



^2,odd\ ^J 2CJ\f\ lf2Md 



(B.29) 



which in terms of V'2,odd alone becomes 

Vi'2,odd = 2\CN'£^M ^2,oAd (A = -z"). (B.30) 

Thus, 

^(A) = det(/-2A£A/'£'^A^). (B.31) 

The previous argument doesn't quite work for dealing with the zeros of A(A), 
since C is singular and one therefore cannot draw the conclusion that the matrix 
LME^ M. appearing in (B.31) is oscillatory (only totally nonnegative, which 
is not enough to show simplicity of the zeros). However, a slightly modified 
argument does the trick. Note that the first row and the last column of the 
2K X 2K matrix in (B.29) are zero. Thus '\\)2(x\) = if (B.29) is satisfied, and 
the value of '^■i{x2K) doesn't really enter into the problem either (it appears only 
in the left-hand side, and is automatically determined by all the other quantities 
in the equation). Therefore (B.29) has nontrivial solutions if and only if there 
are nontrivial solutions to the following truncated [2K — 2) x (2K — 2) problem 
obtained by removing the masses m\ and n^K (i-e., by deleting the first and last 
row and the first and last column): 

f,2,odd\ (B,32) 

V'3,cvcn / 
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where 

■02,0(1(1 = {i'2{x3),-if2{x5), ■ ■ ■ ,^f2{x2K-l)) , 
■03,evon = {i'3{x2) , i^si^A) , ■ ■ ■ , ^f3{x2K -2)) , 

I' = {K -1) X {K -1) identity matrix, 

£' = strictly lower triangular {K — 1) x (K — 1) with £'j = 1 for i > j, 

£' = £ with its first row and last column removed, 
M' =diag(m3,m5,. .. ,m2K^i), 
A/"' = diag(n2, ^4, . . . , 712^-2 )• 

In terms of ■02 odd a-lone, this becomes 

^i,odd - 2A (/' + £')A^'(^')'^-^V2,odd, 
and the conclusion is that 

A{\) = det(/' - 2A (/ + £')-^'(^')^>'V2,odd)> (B.33) 

where (/' + C')M'{E')'^M' is an oscillatory {K - I) x {K - 1) matrix (by the 
previous argument). This shows that A{\) has positive simple zeros too. D 

B.3 Expressions for the coefficients of A and A 

From equations (B.27) and (B.31) we can extract nice and fairly explicit represen- 
tations of the coefficients of the polynomials A{\) and A{\). These coefficients 
are of particular interest, since they turn out to be constants of motion for the 
peakon solutions to the Geng-Xue equation. (It is not hard to show, using 
the Lax pairs, that ^(A) and A{X) are independent of time; the details will be 
published in a separate paper about peakons.) 

First a bit of notation: (j.) will denote the set of fc-element subsets of 
a set 5', and [K] is the set {1,2, ..., iiT}. For a matrix X and index sets 
I = {ii < ■ ■ ■ < im} and J = {ji < • • • < j„}, we write Xjj for the submatrix 
obtained from X by taking elements from the rows indexed by / and the columns 
indexed by J; in other words, Xjj — {Xi^j^)a=i,...,m- 

b=l....\n 

To begin with, (B.27) says that A{X) = det(/ - 2A (/ + C)M£Af), which 
shows that the quantity [A]k from (B.14) (the coefficient of (— 2A)'^ in A{\)) 
equals the sum of the principal k x k minors in (/ + C)A4£J\f: 

[A]k^ J2 det{{I + C)M£M)jj. (B.34) 

A general fact is that for any K x K matrix X and for any fixed J G ( j, ), we 
have the identity 

det{{I + C)X)jj^ Y^ detX/j, (B.35) 

with summation over all index sets / of size k that are "half-strictly interlacing" 
with J: 

14 J ^=4> ii<ji<i2<J2<---<ik<jk- (B.36) 
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(This is similar to, but much simpler than, the "Canada Day Theorem" about 
certain sums of minors of symmetric matrices, which appeared in the context 
of Novikov peakons [17, 15].) Equation (B.35) can be proved by expanding 
det((/ + C)X)jj with the Cauchy-Binet formula and computing the minors of 
/ + £ by applying the Lindstrom-Gessel-Viennot Lemma to the planar network 
for / + £ in Figure 5. (We briefly recall the statement of this lemma: if X is the 
(weighted) path matrix of a planar network G, then the minor det Xjj equals the 
number of vertex-disjoint path families (or the weighted sum over such families) 
connecting the sources indexed by / to the sinks indexed by J.) Alternatively, 
one can do row operations directly, as follows: 



det((/ + C)X)jj = det(((/ + C)X),^,X 



3t 

r,s— 1 



det(^^ X^j 

771—1 

det ( V" X^j^ 

V — / r,s— 1 

h 32 jk 

E E ■•• E det(x„, 

2l = l'2=jl+l «fc=ifc-l+l 



k 



^ det Xi,j 



jl<i2<J2 

i2<i3<33 



jk-l<ik<jk 



(In the second line, we used the definition of C. In the third line, we have 
subtracted from each row the row above it; jq = by definition. Next, the 
summation index is renamed from m to ir in row r; this lets us use multilinearity 
to bring the sums outside of the determinant.) Applying this fact to (B.34), we 
obtain the desired representation 

[A]k= Y. dct{N£M)ij, (B.37) 

This is useful, since these determinants can be evaluated using the Lindstrom- 
Gessel-Viennot Lemma on the planar network shown in Figure 6; see Example B.2 
below. 

In an entirely similar way one derives the identity 

det{CX)jj= J2 detX/j, (B.38) 

with the other type of "half-strictly interlacing" relation 

I^J ■i=> ii< ji<i2< j2<---< Ik <jk- (B.39) 
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Figure 6. A weighted planar network (illustrated in the case K = A) for which 
MEN is the weighted path matrix; this means that the (i,j) entry is the weighted 
sum of all paths from source i to sink j, each path being counted with a weight 
equal to the product of its edge weights. 



Indeed, we can use the network for C in Figure 7, or do row operations: 

k 

r,s— 1 
k 

7',S— 1 

k 
mj. 



det{CX).,j = det[{CX)j^j^ 



Jr-l 

det( Y^ X 

7n— 1 

= detf Y^ X. 



Ji = l J2=ii ik=jk-i 

= Y detX/j. 

n<ii 
ii<i2<i2 

i2<i3<i3 

Jk-l<ik<jk 

From (B.31) we then see that the coefficients defined by (B.22) are given by 
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Figure 7. A planar network (illustrated in the case K = 4) for which £ is the 
path matrix. There is one path from source node i to sink node j ii i > j and 
none otherwise. 



which by the identity above amounts to 

[I]fe= J2 dct{m^M)ij^ Y. det(M£:AA)j,. 



(B.40) 









Again, we can use the network for MEM in Figure 6 to read off these determinants 
(but note that the transposition in the last step of (B.40) has the effect that 
now the sources are index by J and the sinks by /). 

Example B.2. Consider the case K — A. We compute the coefficients [A\k 
using (B.37) and the planar network in Figure 6. 

Note identities such as £'12 £'23-^34 = £14 and {l-El^) + E2z-{l~El2)-E2z = 
1 — _B^3, which are used repeatedly when computing path weights. For example, 
the determinant corresponding to / = {f,3} and J = {f ,4} is found by locating 
all pairs of path connecting source f to sink f and source 3 to sink 4 in Figure 6, 
and having no vertices in common. There is only one path f -> f , and its weight 
is mi ■ 1 ■ E12 ■ n2- Then there are three paths 3 — >■ 4 not touching this first path, 
and their weights are 

ms • (f - Ej^) ■ £56 • Eer ■ £78 • ns, 

m^ ■ E45 ■ (f — £34) • £45 • £56 • Eqj ■ Eys ■ ng, 

ms • £45 • £34 • (f — £23) • £34 • £45 • £56 • £-67 • £78 • "8j 

or, in other words, 

1715(1- El^)E5sns, 

W5(-B45 - £^35)-E'58?^8, 

m5(£|5 - E^5)E5sns. 

Multiplying each of these by the first weight miEi2n2 gives the weights of 
the three vertex-disjoint path pairs 13 -^ 14, which we add up to obtain the 
determinant (according to the Lindstrom-Gessel-Viennot Lemma): 

det(7W£A/')i3,i4 = miEi2n2 ■ m^il - E^^)E5sns. 
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The formulas for the coefficients [A]k found in this way are 

(B.41a) 
= miEi2n2 + miEi^n^ + miEiQUQ + miEisUs + msE^in^ 

+ m-iE^Qna + mgi^agrig + m^E^QUQ + m^E^^jis + mrE7sns, 

il<jl<i2<J2 

= det(M£AA)i2,i2 + det(7W£AA) 12.13 + det(Mf AA)i2,i4 
+ det(7W£:A/')i3,i3 + det(7W£7V)i344 + det(Al£7V)i4,i4 
+ det(7W£:AA) 13,23 + det(7W£A/') 13,24 + det(Al£A/') 14,24 
+ det(M£A/') 14,34 + det(7W£A/') 23.23 + det(X£A/')23,24 
+ det(Al£7V)24,24 + det(Al£7V) 24.34 + det(Al£7V)34,34 

= miEi2n2 ■ 777,3(1 - El^)E34n4 + miEi2n2 ■ m3(l - El^)E'i(inQ 
+ ■miEi2n2 ■ 7773(1 - E23)E3sns + miEi2n2 ■ 7775(1 - El^)E^(in(i 

+ 7771^12772 • 7775(1 - El^)E5^n^ + 7771^12772 • 7777(1 - El.j)E^^ns, 
+ 777li?l4774 • 7775(1 — El^)E5QnQ + 777li?l4774 • 7775(1 — Elr,) E<ig,n^ 

+ miEiifii ■ 7777(1 - El^)EYsns + miEiene ■ 7777(1 - El.j)E7sns 

+ msEsiU^ ■ 7775(1 - El^)E5enQ + m^E-ani ■ 7775(1 - El^)E^s.ns, 

+ mzE^iTii ■ 7777(1 - El^)Ersns + msEsene ■ 7777(1 - EQ.j)E7sns 

+ m^E^QUe ■ 7777(1 - EQ^)Ersn8, 

(B.41b) 

il<jl<i2<J2<i3<J3 

= det{MSAf)i23,i23 + det(7W£7V) 123,124 + det(A^£7V)i24,i24 
+ det(7W£:A/') 124,134 + det(A^£A/') 134,134 + det(A^£A/') 134,234 

+ det(A^£7V)234,234 

= miEi2n2 ■ 7773(1 - £;23)£;34774 • 7775(1 - E%)E5ene 

+ miEi2n2 ■ 7773(1 - E23)E34n4 ■ 7775(1 - El^)E5sns ^ ' '^' 

+ 7771^12772 • 7773(1 - El^)E3ini ■ m'j{l - El.j)E^^ns, 

+ miEi2n2 ■ 7773(1 - E2-i)E3ene ■ 7777(1 - EQ.j)E7sns 
+ miEi2n2 ■ 7775(1 - El^)E^finQ ■ 7777(1 - EQ.j)E7sns 
+ miEiifii ■ 7775(1 - El^)E5ene ■ 7777(1 - EQ.j)E7sns 
+ msEs^Ui ■ 7775(1 - El^)E5ene ■ 7777(1 - EQ.j)E7sns, 



and 



[A]4 = Yl <^(ii{M£J^)^l^2^3^4Jlj2J3J4 

il<jl<i2<J2<i3<J3<ii<ji 

= det(Al£7V)i234,i234 (B.41d) 

= 7771^12772 ' 7773(1 - ^23)^34774 

• 7775(1 - El^)E5Gne ■ 7777(1 - Elj)E'js,ns,. 
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(Cf. the expressions (B.16) for the lowest and highest coefficients [A]i and [A]k 
in general.) 

Similarly, we compute the coefficients [A]k using (B.40): 



[A]: 



[A]i = ^det(Al£AA)j, = ^(Al£AA)j-, 

= m3i?23"-2 + 'rnvjE2r,n2 + TO7i?27"2 
+ m5E45n4 + miE47n4 + mrEQinQ, 

J2 det(Al£A/')j,j,,,,,, 

il<jl<i2<J2 

= det(7Wf A/')23,i2 + det(7W£:AA)24a2 + det(A^£:7V)24,i3 

+ det(>l£A/')34,i3 + det(7W£A/')34,23 
= m3,E2zn2 ■ mzE4^{l - E^4)n4 

+ ™3-E^23"2 ■ rriTEij^l - El^)n4 

+ ™3£^23"-2 ■ mrEer^l - E^Q)ne 

+ ■m5E25n2 ■ mrEej{l ~ Elf^)nQ 

+ mzEi^n^ ■ rriTEQiil - EIq)!!,^ 



and 



[A]3 = ^ det(X£AA)j-^j,j-3,,^,,,3 

il<jl<i2<J2<i3<J3 
= det(A^£A/')234423 

= m3^23"'2 • m^Eizil - El^)n4 ■ mjEej{l - E^Q)nG. 
(Cf. the expressions (B.23) for [A]i and [AJa'-i in general.) 

C Guide to notation 



(B.42a) 



(B.42b) 



(B.42c) 



For the convenience of the reader, here is an index of the notation used in this 
article. 



m{x), n{x), z, *(x;z) = (V^i, -02, V'3)'^ 
Spectral problems for '^ 

g{y), h{y), A = -z^, $(y; A) = (0i, ^2, </>3)^ 

Spectral problems for $ 

Coefficient matrices A{y; A), A{y; A) 

J = ( -1 j 

Involution X{Xy = JX{-X)-^J 

Fundamental matrices U{y; A), U{y; A) 

Transition matrices 

5(A) = t/(l;A), 5(A) = [7(1;A) 
Weyl functions 

M^ = -52i/53i, ^=-5n/^3i 



Section 1 


(1.3a) + (2.4) 


(1.4a) + (2.4) 


(2.1) 


(2.2), (2.3) 


(2.6) 


(2.8) 


(2.9) 


(2.10), (2.11) 


(2.12) 


(2.17) 
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< XN 



+ mN-iSxr. 
- riNSxj^ 
N ^2K 



ha 
Vk 



2n2a cosh X2a 



Twin Weyl functions 

W = -S21/S31, Z = -Sn/Ss, 
Bilinear form (<i?,ri) = J_^ ^{y)'^ J n{y) dy 
Adjoint spectral problems for 

n = n{y;X) = {uJi,uj2,^3) 
Adjoint Weyl functions 

' W* = -S32/S31, Z* = -533/^31 
Twin adjoint Weyl functions 

W* — —S32/S31, Z* — —S33/S31 
Discrete interlacing measures 

m = miSx^ +11136x3 

n = n2Sx2 + n^Sx^ + 

with Xi < X2 < 
Transformed measures 

9 = 9i^vi +925y3 H 1 

h = hi5y^ + /i2(5y^ H 

with i/fc = tanh Xk , 

ga = 2TO2a-lCOsha;2a-l, 

Interval lengths l^ — yk+i 
Propagation matrices 

r / A N / 1 \ r^" 
Lk(X} = 10, 

Transition matrix in the discrete case 

S{\) = L2K{\)[\^]L2K-liX)[°J--- 

and its partial products Tj{X) 
Twin transition matrix 

SiX) = L2K{X)[l]L2K-liX)['^]--- 

and its partial products Tj{X) 
Eigenvalues 

= Ao < Ai < •• • < Ak 

= /io < Ml < • • • < fJ'K-l 

Residues of Weyl functions 

a^, bj, 600, c,, dj {l<i<K,l<j <K-l) 
Spectral measures 

" = E^=i "i'^A, , /3 = Y.f=i^ bjSf, 
Weyl functions as integrals 



1 X 2^y 
1 y 
1 



Jfij 



da(x) 



Entries of r(A) = T,(A) (for some fixed j) 

Q ^ ~T32, P = T22, R = ^2 
{W w P/Q, Z w R/Q) 
Residues of adjoint Weyl functions 

Adjoint transition matrix 

5*(A) = 5(-A)-i = J5(A)^J 
Adjoint Weyl functions in terms of S* 

w* = +s*2jsi^, Z* = -^1*1/^3*1 

Moments and bimoments of spectral measures 
ak ^ j x^da{x), l3k = J y''dl3{y), 



(2.17) 
(2.20) 

(2.22), (2.25) 

(2.24) 

(2.24) 

(3.1) 

(3.3) 
(3.2) 
(3.4) 
(3.5) 

(3.7), (3.8) 



(3.9) 




(3.10) 




(3.23) 




(3.24) 




Theorem 


3.8 


Theorem 


3.10 


(3.43) 




(3.44) 




(3.45) 




(3.52) 




Theorem 3.18 


Remark 


3.19 


(4.1); see 


also (A.l) 
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Determinant /C„ involving binioments 
Spectral map 

Pure peakon sector V C R*^ 

Admissible spectral data TZ C R"*^ 

Forward map S : V ^ TZ 

Inverse map T: 7?. — > P 
Cauchy biorthogonal polynomials Pn{x), qn{y) 
Vandermonde-type expression 

T{x) =r(a;i,...,x„) = Y{^^J{x^ + XJ) 
T{x;y) ^T{xi,...,Xn;yi,.--,ym) 

= n:unr=i(-^+y.) 

Generalized Heine-type integrals 

^rs _ r A(a:)^A(y)^(];]='^.)''(ny^)°'^""(^)'^^"'fa) 

'Jnm J T{x;y) 

where cr„ = {x & R" : < xi < • • • < a;„} 
Degenerate cases JJ'^, J^o- 'To'o' 
The basic bimoment determinant 

Dn = det(iij)jj-^g — j„„ 
General discrete setup 

a = Ell a^^A, , /3 = Ef=i ^i'^M, 
(v4 = X, S = A' — 1 in the main text) 
Heine-type integrals as sums in the discrete case 

JnL = E/6([;Ji) E,/e(i^i) "^IJ ^"i^-i t^J^J 
where 

[^] = {1,2,...,^} 

in) ~ ^^^ ^^ n-element subsets of [A] 

X'jai fi'jbj = (riie/ ^ra») (rije,/ A^l^jJ 



* 



A?A2, 



/,/ 









T/j = r(Ai^, . . . , Aj„;^j 
A/i/2 = lliie/i, 126/2(^^1 - ^^2) 
J^nm written out in the case A = 3, 

Polynomials Ak{X), Sfe(A), Cfe(A) 

Jump matrix Sk{X) 

(A(A),S(A),C(A)) = (A^(A),S^(A),Cjv(A)) 

Coefficients [A]„ [S]„ [C], in A, B, C 

Eab = e-l^»--jl 

Polynomials A^(A), Sfe(A), Cfe(A) 

Jump matrix 5fc(A) 

{A{X),B{X),d{X)) = {AN{X),BN{X)^CN{X)) 
Coefficients [A]„ [B]„ [C], in A, B, C 

1p2,cvcn, '03,odd, C, £, M, N 



(4.2); see also (A.2) 
(4.5); see also (A.26) 



Definition 4.7 

(A.4), (A.5) 



(A.12) 
(A.13) 

(A.14) 



(A.15) 
(A.16) 

(A.17) 

Section A.4 

(A.31) 



(A.33) 




(A.34) 




(A.35) 




Examp 


e A.2 


Lemma 


A.3 


(B.8), (B.IO) 


(B.ll) 




(B.12) 




(B.14) 




(B.15) 




(B.17), 


(B.18 


(B.19) 




(B.20) 




(B.22) 




(B.25) 
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"02,0(1(1, "03, oven 

V-^.ovon, ^Udd, -C', £', M', N' 

[K]^{l,2,...,K} 

(^) , the set of fc-element subsets of a set S 

Index sets I = {ii < ■ ■ ■ < im\, J = {ji < • • • < jn\ 



Submatrix Xjj = [Xi^j^)a=i,...,m 

b=l,....,n 

"Half-strictly interlacing" relations: 

I ^ J ^=> ii < ji < 12 < J2 < 
I A J ^=^ ii < ji <i2< J2 < 



<ik< jk 
<ik < jk 



(B.29) 
(B.32) 



Section B.3 



(B.36) 
(B.39) 
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